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CHAPTER 1 

GENERAL INFORMATION 



THE FINITE ELEMENT CODE AT I LA HAS BEEN DESIGNED BY THE ACOUSTICS 
LABORATORY OF THE "INSTITUT SUPERIEUR D' ELECTRON I QUE DU NORD" 
(I SEN. LILLE. FRANCE) FOR THE "GROUPE D* ETUDES ET DE RECHERCHE EN 
DETECTION SOUS-MARINE" (GERDSM. DCAN. FRENCH NAVY. TOULON). THE 
SOFTWARE ENGINEERING IS PROVIDED BY THE SINAPTEC COMPANY (LILLE. 
FRANCE). 
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1.1 INTRODUCTION 



The finite element code AT I LA has been specifically developed 
to aid the design of piezoelectric devices, mainly for sonar 
applications. Thus, It Is able to perform the modal analysis of 
both ax Isymmetr lea I and fully three-dimensional piezoelectric 
transducers. It can also provide their harmonic response under 
radiating conditions: nearfleld and farfleld pressure, 
transmitting voltage response, directivity pattern, electrical 
Impedance, as well as displacement field, nodal plane positions, 
stress field and various stress criteria... Its accuracy and Its 
ability to describe the physical behavior of various transducers 
(Tonpllz transducers, double headmass symmetrical length 
expanders, free flooded rings, f lextens Iona I transducers, bender 
bars, cylindrical and trilaminar hydrophones...) have been checked 
by modelling more than twenty different structures and comparing 
numerical and experimental results. These tests were~'all 
successful, provided that sufficient care was previously given to 
the physical meanings of the models. Most of them, as well as the 
general modelling rules, have been described in reports, papers 
and communications, which are partially listed In the next 
sect ion. 

AT I LA has been designed for the " Groupe d* Etudes et de 
Recherche en Detection Sous-Marine, G.E.R.D.S.M." from the 
"Direction des Constructions et Armes Navales, D.C.A.N.", In 
TOULON (FRANCE), by the Acoustics Laboratory of the "Inst I tut 
Superleur d'ElectronIque du Nord, I.S.E.N.", In LILLE (FRANCE). It 
has first been a research tool, for specific optimisation 
purposes, but, due to the request for Its use by a large number of 
people. It has then been made more user friendly and documented, 
while a specific effort has been carried out to Improve Its 
versatility as well as to provide more efficient graphic display 
abilities. At the present time, ATILA has been run on various IBM, 
VAX and UNiVAC mainframes, as well as . a FPS array processor and 
VAX and MASSCOMP workstations. Nevertheless, ATILA Is not a "black 
box" and users are assumed to have a firm background In fJnIte 
element analysis. 

ATILA has been developed over the past 8 years by a lot of 
people, but mainly by K. ANIFRANi, R. BOSSUT, J.L. CARTON, J.C. 
DEBUS, J.N. DECARPIGNY, B. DUBUS, R. EDDE, B. HAMONIC, D. MOREL, 
P. TIERCE at I.S.E.N., and B. TOCQUET and D. BOUCHER at 
G.E.R.D.S.M.. Current developments. Including new algorithms, new 
material description, material loss description and finite 
element/Integrai equations coupling are proceeding at I.S.E.N., 
while software engineering Is now provided by the SINAPTEC 
company, LILLE (FRANCE). 
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Information and help about AT I LA can be obtained, under 
specific conditions, from: 

. J.N. DECARP IGNY 

Acoustics Laboratory 

Inst I tut Super leur d'ElectronIque du Nord 
41 boulevard Vauban 
59046 LILLE Cedex. FRANCE 

Tel: (33) .20.30.62.20 TeleX: GEFIRN 120.369 

. P. TIERCE 
SINAPTEC Company 
. 41 boulevard Vauban 
59800 LILLE. FRANCE 

Tel: (33) .20.57.80.04 Telex: GEFIRN 120.369 



Following this section and the reference list, the reader will 
find the description of the successive stages of an ATILA Job and 
their connection to the manual organisation. Then, detailed 
Information Is collected In five chapters. 
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1.3 GENERAL ORGANISATION OF AN AT I LA JOB 



Running an ATILA Job requires different well defined stages. 
This section describes briefly their organisation and their 
connections to the various following chapters of this manual. 

1.3.1. THE MODEL DEFINITION. 

Being given the physical problem, the user has to define the 
type of analysis which Is required: computation of the static 
deformation of an elastic stucture under lumped mechanical load, 
modal analysis of a piezoelectric stack with short circuit 
electrical condition, computation of the transmitting voltage 
response of a transducer... This type of analysis has to be 
related to solving algorithms which are listed and described In 
chapter 2. New solvers are currently under test or development In 
I.S.E.N. Acoustics laboratory and detailed Information c¥n' be 
provided upon request. The user has also to determine the type(s) 
of elements which Is (are) required: three-dimensional elements, 
ax I symmetrica I , plane stress or plane strain elements, plate or 
shell elements relying on specific approximations... The ATILA 
element library Is described In chapter 4 . 

At this stage, a first estimation of the computation size can 
generally be made, using a gross determination of the node and 
degree of freedom numbers and previous CPU time measurements or 
classical numerical rules. This estimation Is Important since It 
can avoid the generation of unrealistic meshes. At this stage 
also, the user has to gather the requested physical parameters, as 
described In chapter 4, and to ensure their accuracy, which very 
definitely determines the accuracy of the result. Though not 
directly Involved In the data file generation, this stage Is 
fundamenta I . 

1.3.2. THE MESH GENERATION. 

-Following the previous choice of element type(s), as de-fined 
during the first stage, the mesh generation Implies the splitting 
of the whole domain Into a given number of elements. This Is 
realised by defining nodes and assigning them to corresponding 
elements. Node definition Is made by listing the node coordinates 
In a given order, called the node numbering order. Then, 
assignment of the nodes to the corresponding elements Is obtained 
by listing for each element. In a given order named Its topology, 
the nodes which belong to It. The node coordinate description Is 
detailed In chapter 3, together with all other entries. The 
topology description Is presented also In chapter 3, after the 
node coordinate description, but more details are available on an 
element by element basis In chapter 4. 
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During the mesh generation, the user can use the Isoparametric 
character of all the ATILA elements, which allows easy curved line 
modelling and lower density meshes. Nevertheless, a great deal of 
attention has to be paid to the element geometrical aspect, as 
well as to the mesh step which has to be related, for example, to 
the minimal wavelength value of Interest. 

Mesh generation can be fully carried out by the user. However, 
in many cases, partial regularity of the mesh allows the use of an 
automatic generation which provides node coordinates and element 
topologies. Then, the user has only to define a gross splitting, 
built with super elements, and to select their automatic splitting 
Into finite elements. This easier procedure, described in chapter 
6, relies upon the use of the MOSAIQUE code. 

1.3.3. THE DATA FILE PREPARATION. 

An ATILA data file includes entries, specific data and 
comments, it Is devoted to the description of the type of 
analysis, the node coordinates, the element topology, the material 
properties, the element geometrical properties, which are all 
entries, and the loading and boundary conditions which are 
formatted data, in some cases, it also Includes specific mixed 
finite eiement/plane wave model data and graphic display data~. The 
whole data file preparation Is described In chapter 3 which has to 
be carefully read. Specific element information which are required 
during this preparation can also be found In chapter 4. Finally, 
when the automatic mesh generation code MOSAiQUE Is used, data 
file preparation can be, in a large part, automatically obtained, 
as explained In chapter 6. 

A detailed check of the data file can be processed by running a 
specific progam, named PGRAPH, which produces a complete listing 
of the various data. Including the generated degree of freedom 
numbering which reflects exactly the prescribed boundary 
conditions, as well as different kinds of graphic displays (total 
or partial mesh, with or without element, node or degree of 
freedom numberings). The use of PGRAPH Is described In chapter 5. 

1.3.4. THE RUNNING OF THE CODE. 

Running of the code requires the edition of a main program 
which has to set the array sizes, call successive subroutines to 
read the data file, compute the elementary stiffness and-, mass 
matrices, assemble these matrices Into the global stiffness and 
mass matrices, solve the equations and display the results. This 
progam generation is automatic, given the corresponding data file, 
and is provided by running the program PGEN. PGEN will Issue 
successive questions, on the terminal screen, to specify various 
user's requests, such as array storages for a delayed solving, 
array storages for a postprocessing, graphic displays... it will 
provide various files. Including the main program FORTRAN source 
and, in the case of a VAX computer under VMS operating system, a 
COMMAND file to compile this main program and edit the links. 
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Running of the code will provide a detailed result listing and, 
according to the specific requests, graphic displays or stored 
arrays for restarting another solving or for postprocessing. It 
has to be noted that PGEN Is also able to generate specific main 
programs to restart given runs. Every Information which concerns 
the program running as well as PGEN Is contained In chapter 5. 

1.3.5. RESULT USING AND POSTPROCESSING. 

As soon as the run Is completed, a LISTING file Is available 
which contains all the results: elgenfrequencles and eigenvectors, 
displacement and pressure fields, electrical Impedance... If 
requested, specific arrays are also stored for postprocessing, and 
mainly the array PLO which contains displacement and pressure 
fields for every loading case, mode number or excitation 
frequency. Graphic displays of the distorted structure can be 
easily obtained on various workstations using PGRAPH . • PGRAPH will 
Issue FORTRAN and COMMAND files which rely upon the GKS package 
and have to be executed to provide the display. As mentioned 
previously, the use of PGRAPH Is described In chapter 5. 

1.3.$. SUMMARY. 

As a summary, the next page displays a simple flow chart “of an 
ATILA finite element modelling. A large part of this flow chart 
will be reproduced in a more detailed and accurate fashion In 
chapter 5, with reference to the various files and programs. 
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2.1 GENERAL FORMULATION 



This section describes the genera i formuiation which has been 
retained to model radiating piezoelectric transducers using the 
ATiLA code and provides a complete list of the available types of 
analysis. Then, for each type of analysis, the following sections 
give more information and specific details. 




Boundary 
surface S. 



To model a radiating piezoelectric transducer using the ATiLA 
code, the finite element mesh must include the structure as. well 
as a part of the fluid domain, as described on the above figure. 
Thus, the unknown fields are the displacement field u in the whole 
structure, the electrical potential v in the piezoelectric 
material and the pressure p in the fluid. Then, the equations 
which are solved are, first, the equation of motion in the 
structure, second, the Poisson equation in the piezoelectric 
material and, third, the Heimhoitz equation in the fluid. 
Kinematic and dynamic continuity equations between displacement 
and pressure fields are prescribed at the interface, due to the 
variational formuiation, and an appropriate damping condition Is 
applied to the spherical external surface S that surrounds the 
fluid doma I n . 
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In matrix form, the complete set of equations associated with 
this problem Is: 



![Kuu] - w2[M] 
1 


I— 1 

7s 

< 


- CL] I 

1 


yi 

•1 


F 


! [Kuv]"^ 

1 


[Kyy] 


[O]''' ! 

1 


1 

VI - 
1 


-Q 


! -a2c2w2[L]''’ 


[0] 


[H] - w2[Mi]| 


1 

PI 


ac^F 


where : 










- U, F, V, Q 


and P are 


vectors that 


contain 


the 



_ _ _ _ _ nodal values of 

the displacement field, the lumped applied forces, the 
electrical potential and charges, the pressure field, 

F' Is a vector that contains the nodal values of the outgoing 
flux through S, 

[Kuu] and [H] are the solid and fluid stiffness matrices, 

[M] and [M-|] are the solid and fluid mass matrices, 

[Kq^] and [Kyy] result from the assembling of the piezoelectric 
and dielectric matrices, 

[L] Is the fluid-structure Interface connectivity matrix, 

[0] Is a null matr lx, 

a and c are the fluid density and sound speed, 
w Is the circular frequency, 

T means transposed. 



In fact, 
assoc I ated 
Moreover , 
electr leal 
fol lowing 
potential 
preceeding 
AT I LA code 



Internal electrical potential degrees of freedom are 
with null nodal charges and can be easily condensed, 
nodes which belong to hot electrodes have the same 
potential degree of freedom which can be factored, 
simple algebraic operations. Finally, the reference 
can be arbitrarily put equal to zero. Thus, the 
matrix equation Is automatically transformed by the 
Into the following form: 



[Kuu] - w^[M] 


!S uv 


- [L] I !U1 IF 

III 1 


uv^ 


> 

> 


0 1 |V| - l-l/Jw 

III 1 


-a^c^w^CL]]’ 


0 


[H] - w2[Mi]| !p! lac^F* 



where elastic, piezoelectric and dielectric terms are modified by 
the Internal electric potential static condensation, K ^Y and 0 
are now vectors, Kyy Is now a scalar, V is the applied electrical 
potential and I Is the electrical current entering Into the 
structure. Then, If a dipolar damping condition Is applied on the 
external boundary surface S, F' Is given by: 



F* - -(1/ac).(1/R + Jw/c).[D] P 

+(1/ac) .{(1/R - Jw/c)/(1 + w2r2/c2)).[D'] P 
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where R Is the radius of S, [D] and [O'] are obtained by 
assembling radiating surface finite elements from S, the first 
term is the monopolar contribution associated with a spherical 
wave Impedance and the second term is the dipolar contribution. 

Starting with this general problem, a lot of different specific 
analysis can be processed. Those for which AT I LA provides a solver 
are the fol lowing: 

-computation of the displacement field of an elastic structure 
under a static lumped or distributed load (STA 1), 

-computation of the displacement field and the electrical 
potential for a piezoelectric structure under a static lumped or 
distributed load (STA 2), 

-modal analysis of an elastic structure, providing the 
eigenfrequencles and eigenmodes (MOD 1), 

-modal analysis of a piezoelectric structure, providing the 
eigenfrequencles and eigenmodes under resonance or ant I resonance 
electrical conditions (MOD 2), 

- modal analysis of a closed fluid domain, with zero pressure 
and/or zero flux boundary conditions, providing the 
eigenfrequencles and eigenmodes ( MOD 3), 

-modal analysis of a closed hydroeiastic system, built with two 
coupled elastic and closed fluid domains, providing the 

eigenfrequencles and eigenmodes (MOD 4 ), 

-modal analysis of an elastic structure flooded In an 

Incompressible fluid (added mass effect only), providing the 
eigenfrequencles and eigenmodes (MOD 5), 

-harmonic analysis of an in-aIr piezoelectric structure, 
providing its electrical Impedance and displacement field for a 
given excitation frequency (HAR 1), 

- harmonic analysis of an In-air piezoelectric structure, using 
the mixed plane wave/fInIte element method and providing the 
displacement field for a given excitation frequency (HAR 2), 

- radiated pressure computation for a given displacement field of 
a vibrating structure at a given frequency (HAR 3), 

-harmonic analysis of a radiating piezoelectric transducer, 
providing the electrical impedance and the displacement and 
pressure fields in the whole mesh for each given excitation 
frequency (HAR 4 ), 

- harmonic analysis of a radiating piezoelectric transducer, using 
the mixed plane wave/finite element method and providing the 
displacement and pressure fields In the whole mesh for each 
given excitation frequency (HAR 5). 
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2.2 STATIC ANALYSIS OF AN ELASTIC STRUCTURE (STA 1) 



In this case, only the displacement field Is concerned, while w 
has to be put equal to zero, and thus the whole matrix equation Is 
reduced to: 



[Kuu] y - F 



The user has to provide the loading F and the solver computes jJ. F 
has to be defined strictly for the finite element model, which 
means that, for example. It has to be divided by 2*pl If the model 
is ax Isymmetr leal , by 2 If It Is applied In a symmetry plane 
limiting the mesh... A frontal assembling Is performed, coupled 
with a static Gauss condensation, and then a classical 
backsubst I tut ion algorithm restores the whole response. Several 
different loading cases can be solved simultaneously. The user has 
to take care of the presence of rigid body modes, which are 
associated with singularities of the system of equations. FlhaTly, 
If distibuted loads are applied, the user has to compute the 
corresponding lumped forces, using the element Interpolation 
funct Ions. ' 
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2.3 STATIC ANALYSIS OF A PIEZOELECTRIC STRUCTURE (STA 2) 



In this case, the matrix equation is reduced to: 



I [Kuu3 


K UV' 

1 


!U| 
1 1 


1 FI 


! K T 

1 N UV 


Kvv ! 


1 1 

1 V 1 
1 V 1 


!-Q! 



If the external electrical connection Is a short circuit, V is 
equai to zero. The proceeding equation Is reduced again to: 



[Kuu] U “ f 



The user has to provide the loading F and the solver computes U, 
following the same procedure as In the proceeding section. If the 
external electrical connection Is open, Q Is equal to zero an^:' 



CKuu3 5 uv! Sy! !£! 

I I I ^ I I 

K uv^ Kvv ! !vl !ol 



The user has to provide F and the solver computes U and the 
induced eiectricai potent ia*l V. Here again, attention has to be 
paid to rigid body modes as well as. In the second case, to the 
uniqueness of the eiecticai potential. As for STA 1, several 
loading cases can be processed at the same time. 
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2.4 MODAL ANALYSIS OF AN ELASTIC STRUCTURE (MOD 1) 



In this case, only the displacement field Is of concern and the 
loading force Is put equal to zero. The matrix equation Is reduced 
to: 



CCKuu] - w 2[M]) y - 0 



Eigenvalues and eigenvectors from this linear system of equations 
are the el gen frequencies and elgenmodes of the structure. 

The AT I LA code uses frontal assembling coupled to a Guyan 
condensation. Thus, the user has to provide the number of master 
degrees of freedom which have to be retained In the final 
diagonal I sat ion phase. To do this, he has to keep in mind that the 
CPU time which Is reqested for the d lagona I I sat Ion phase Is 
roughly proportional to the cube of the master degree number, 
absolute value of this time depending upon the algorithm, the 
required accuracy, the computer and the compiler option. He has 
also to keep In mind that accurate e Igenf requenc les are obtained 
only for: 



w2 « Min(K| |/M| I ) 



where K| | and M| | are the diagonal terms of the stiffness and mass 
matrices associated with the condensed slave degrees of freedom. 
Thus, he has to perform a compromise between these opposing 
requirements. For a given master degree number, the choice of the 
specific degrees which have to be retained can be provided by the 
user or automatically selected by the code, using the HenscheM's 
cr i ter i a . 

The diagonal I sat ion algorithm which Is used by ATILA relies 
upon the classical Householder tr I d I agona I isat Ion method, coupled 
to a Sturm sequence, it can be processed in simple or double 
precision. CPU times, as obtained on an IBM 4341 Mod. 2 (4-.. MO), 
are reported on the graphic displays of pages 8 and 9. Recently, a 
new d I agona I isat Ion algorithm has been Implemented. It provides 
the e I genf requenc les In double precision and the eigenvectors In 
single precision. Though much more efficient, it Is still under 
test and can be only used by advanced users, following a short 
modification of the ATILA main program*. In ail the cases, ATILA 
provides ail the el genf requenc les, the number of which is equal to 
the number of master degrees, and only the first twenty 
elgenmodes, if available, it provides also the rigid body modes, 
which belong to the previously mentioned first twenty modes. Upon 
the user's request, a backsubst Itut Ion routine can provide the 
slave degrees of freedom for any of the obtained elgenmodes. 
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In some cases, care has to be taken because of mixing between 
modes that have frequencies which are too close, as well as of 
matrix 1 1 l-cond 1 1 Ion Ing. Successive diagonal I sat Ion with different 
precision or different master degree numbers can Illuminate these 
two points. 



*. To use the new d I agona I Isat Ion algorithm, the user has to 
edit the main program FORTRAN source, as obtained for a double 
precision eigenvalue computation, replace CALL D52200 (DLIST,S, 
S,CPDDC,PLO, IP,DM,Z) by CALL S52200 (DL IST,S,CPDDC,PLO, IP ,Z) and 
delete the DM array from the DIMENSION set of Instructions. 
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CPU TIME FOR THE D I AGONAL I SAT I ON PROCESS WITH RESPECT TO THE 
NUMBER OF MASTER DEGREES OF FREEDOM. The parameter is the total 

number of degrees of freedom. 
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2.5 MODAL ANALYSIS OF A PIEZOELECTRIC STRUCTURE (MOD 2) 



in this case, the matrix equation is reduced to: 



I [Kuu] - w2[M] 


K Qvl 
1 
1 


!Ui 
1 1 
1 1 


! 0 1 

m 1 1 


1 T 

1 t uv 


> 

> 


|V| 


l-l/Jw! 



If a resonance condition occurs, V Is equal to zero and: 
CCKuu] - w 2[M]) y - 0 

If an ant I resonance condition occurs, I Is equal to zero and: 
CCKuu] - w 2[M] ” !S uv!S uv ■T/Kvv) y - 0 



In the two cases, all the Information given In the proceeding 
section with respect to the modal analysis of an elastic structure 
is applicable. The only point which has to be added Is that, due 
to the automatic static condensation of the electrical potential 
degrees of freedom, all these degrees of freedom are slaves. 
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2.6 MODAL ANALYSIS OF A CLOSED FLUID DOMAIN (MOD 3) 



In this case, only the pressure field Is concerned. Moreover, 
as zero pressure and/or zero flux boundary conditions have to be 
applied, the F' components which are able to appear In the 
equation are all zero. Thus: 



C[H] - w2[Mii} P - 0 

All the Information provided In section 2.5 with respect to the 
modal analysis of an elastic structure Is applicable. Moreover, 
the user has to remember that the zero flux condition is the 
natural boundary condition of this finite element formulation and 
automatically applies to any boundary. If not superseded by 
another one. This modal analysis can be used to determine the 
Irregular frequencies associated with an exterior Helmholtz 
integral equation. 
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2.7 MODAL ANALYSIS OF AN HYDROELASTIC SYSTEM (MOD 4) 



In this case, structure and fluid domains exist together and 
Interact along their Interface. However, the fluid domain has to 
be closed, no radiation damping being Included In the model. Thus, 
fluid surfaces which are not in contact with the solid structure 
have to be zero -pressure or zero flux surfaces, as Is the case, 
for example, for fluid in an open tank. Restricting the general 
matrix equation of section 2.1 to the displacement and pressure 
field, the equation which is obtained Is: 



[Kuu3 - w2[M] - [L] 

-a^c^w^CL]'*' [H] - w2[M-|] 



iU! 

I I 
I I 



I p I 
I r I 



0 

0 



Apart for a slight difference regarding the algorithm, ail the 
Information given In section 2.4 with respect to the modal 
analysis of an elastic structure Is applicable. With respect to 
the algorithm, the above matrix has to be made symmetrical before 
the ATILA classical d lagona I Isat Ion Is' performed. This 
symmetr isat Ion Is processed using a standard method due to irons. 
However, It has to be pointed out that the symmetr Isat ion method 
requires the inversion of [H] and then implies that no fluid 
Isopressure mode, equivalent to structure rigid body modes, 
exists. Thus, at least a small part of the fluid surface has to be 
a zero pressure surface. Recent developments have led to the 
design of a new solver which avoids this dlfflculty^^ but It Is 
still under test and can only be run by advanced users . Moreover, 
during the assembling phase, as well as for an elastic structure, 
Guyan condensation is performed for both solid and fluid degrees 
of freedom. Care has to be given to the interface degrees of 
freedom which have to be kept in sufficient number. Finally, 
piezoelectric parts can be included in the structure. Then, the 
user has also to provide the corresponding electrical boundary 
cond I t Ions . 



*. To use this new algorithm, the user has to edit the main 
program FORTRAN source, as obtained for this analysis, delete CALL 
R52300 and replace CALL R52200 by CALL S52300 with the same 
parameters. The new algorithm detects automatically the 
isopressure or rigid body modes and deals with them using a 
deflation method, in opposition to the standard algorithm, it can 
also compute the eigenmodes for mere fluid or solid domains, 
without interaction. 
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2.8 MODAL ANALYSIS OF AN ELASTIC STRUCTURE FLOODED IN AN 
INCOMPRESSIBLE FLUID (MOD 5) 



Incompressible fluid condition Is an approximation that Is 
useful to compute the fluid mass loading effect on the low 
frequency modes of an elastic structure. Obviously, It discards 
any radiation effect. If It Is applied to the matrix equation of 
the previous section, this equation becomes: 



CCKuu] - w 2([M] + a2c2[L][H]“‘'[L]''')} U - 0 

where a^c^[L] [H]”^ [L]^ Is the added mass. Then, the classical 
ATILA diagonal I sat Ion provides the structure elgenfrequencles and 
elgenmodes. To obtain this particular condition using ATILA 
requires that all the fluid degrees of freedom are condensed, this 
approximation producing exactly the same result as the 
incompressibility hypothesis . If the fluid domain extension Is 
Infinite, it has to be simply truncated before the mesh 
generation. Successive truncations at various distances from the 
wet structure are the best convergence test. 



After the Guyan condensation of all the fluid degrees of 
freedom, the matrix which has to be diagonal I sed Is already 
symmetrical but the automatic program generator will provide the 
call for the symmetr Isat Ion subroutine because the data file 
corresponds to a fluid-structure modal analysis. Thus, before 
running ATILA the user must edit the FORTRAN file and delete the 
I Ine CALL R52300. 
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2.9 HARMONIC ANALYSIS OF AN IN-AIR PIEZOELECTRIC STRUCTURE (HAR 1) 



In this case, an electrical potential of given frequency Is 
applied to an In-air free piezoelectric structure. Then, the 
general matrix equation of section 2.1 can be reduced to: 



![Kuu3 - 


K Qyl 
1 


!U1 
■ ■ 


1 0 1 
K * ' 


! K T 

1 t uv 


> 

> 


1 1 
IV! 


i-i/jw! 



To obtain this whole equation using the ATILA code implies that 
the applied electrical potential degree of freedom is not 
automatically condensed, which requires a specific Identification 
using the boundary condition data. Then, the solver computes the 
displacement field U, using the first line of the above matrix 
equat ion: 



([Kuu] - w2[M]) y - -K uv.V 



as well as the electrical Impedance Z, using the second line: 
1/Z - -jw (Kvv + K uv''’.y/V) 



The solving algorithm is a simple Gauss elimination, in the same 
run, solving can be performed for 20 different excitation 
frequencies. Moreover, Job restart is possible for new frequency 
values without performing again the whole assembling. The 
Impedance value which Is provided corresponds strictly to the 
finite element model. Thus, the user has to divided the result by 
2’*‘pl If the model is ax Isymmetr lea 1 , 2 If the mesh is restricted 
to half of the structure by a symmetry plane, etc. 
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2.10 HARMONIC ANALYSIS OF AN IN-AIR PIEZOELECTRIC STRUCTURE USING 

THE MIXED METHOD (HAR 2) 



In this case, the driving part of the structure Is modelled 
using the plane wave approximation and a transfer matrix method, 
while the remaining Inactive part of this structure is modelled 
using finite element. These two parts are In contact along a plane 
interface, which exhibits also the plane wave type displacement 
field. Thus, all the degrees of freedom which belong to the finite 
element side of the interface and correspond to a mechanical 
displacement normal to It have to be Identical and are denoted 
hereafter as Uq. Then, the transfer matrix method provides, for 
each given frequency, a mechanical force Fq which is generated by 
the piezoelectric effect and is applied to the finite element part 
at the interface. This force can be written: 



Fq " Zq . Uq + Fy 



Where Zq is a mechanical impedance and F^^ contains the electrical 
excitation effect. Finally, the general matrix equation becomes: 

! [Kuu] - w2[M] 

I Ko'' 

where [Km^] ahd [M] correspond to ail the degrees of freedom apart 
from Uq. The user has to select the excitation frequencies and to 
specify, using the correct boundary condition data, the degree of 
freedom Uq. The solving algorithm is a simple Gauss elimination, 
and the code provides the displacement field in the finite element 
part of the structure. Twenty different excitation frequencies can 
be provided for the same run. Moreover, as in section 2.9, job 
restart Is possible with new frequency values. This mixed method 
is generally well suited tp the modelling of Tonpliz type 
transducers and allows large savings of CPU time . 
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2.11 RADIATED PRESSURE COMPUTATION FOR A GIVEN DISPLACEMENT FIELD 

(HAR 3) 



in this case, a vibrating structure radiates an acoustic wave 
In an Infinite medium and the vibration displacement Is assumed to 
be known for each given frequency. The general matrix equation 
from section 2.1 can be reduced to: 



C[H] - w 2[M-|]) P - ac^F* + a^c^w^ELl^.y 



Then, If monopolar or dipolar damping Is used, following the third 
equation of section 2.1, the damping term can be written under the 
form: 



ac^F* - [G] P 



where [G] stands for a linear complex frequency dependent 
operator. Thus: 



{[H] - w2[Mi] - [G]) P - a2c2w2[L]''‘.y 



The user has to define U and the given excitation frequencies. 
To define U, he has to write a small subroutine, named FABU*, 
which contains the displacement value for each Interface solid 
degree of freedom, and has to be added to the main program. Only 
one displacement field can be defined at a time and this field has 
to be real (phase reference). Ten successive frequencies can be 
processed In the same run and, as for the previous section, 
restart Is available. The solver splits the previous complex 
equation Into two real equations. Then, a matrix inversion 
followed by a Gauss elimination Is performed. For each frequency, 
ATILA provides the real and Imaginary parts of the pressure field, 
printed in two successive tables. 

To obtain accurate results, the user has to give a great deal 
of attention to the damping element choice, if monopolar damping 
Is used, the external boundary surface has to be centered at the 
acoustic center and It has to be In the farfleld. If dipolar 
damping Is used, the external boundary surface has also to be 
centered at the acoustic center but, generally. It can be deeply 
In the nearfleld. However, in this second case, farfleld 
Information Is only available after an appropriate postprocessing. 
Dipolar elements can only be used for ax Isymmetr lea 1 models. 
Three-dimensional dipolar element design Is In progress. 
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Degrees of freedom which belong to the exterior boundary 
surface are automatically master degrees. Those which belong to 
the Interface can be slave degrees but the user has to keep a 
sufficient number of them as master degrees. Finally, care has to 
be taken to keep the maximal distance between two successive fluid 
master nodes In the propagation direction less than one quarter of 
a wavelength. 



The FORTRAN subroutine which has to be added at the end of 
the main program, to provide the given displacement field U, Is 
FABU (CPDDC,U, IDLE , IP) . It Is called by the solving subroutine 
R52400 during the run. IDLE and IP are control Integers which are 
defined In R52400, as well as the CPDDC array. In return, FABU 
provides all the components of U, numbered In the ascending order. 
U has to contain IDT components, where IDT Is the total number of 
degrees of freedom, but. In fact, only the mechanical degrees of 
freedom are concerned. The following lines are provided -..as a 
standard example: 

SUBROUTINE FABU(CPDDC,U, lOLE , IP) 

DIMENSION CPDDCCIP,*), U(») 

U(1) « ... 

U(2) «... 



RETURN 

END 
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2.12 HARMONIC ANALYSIS OF A PIEZOELECTRIC RADIATING STRUCTURE 

(HAR 4) 



This is the general case. Using the operator [G] defined In the 
previous section 2.11, the matrix equation Is: 



ICKuu] - w2[M] 


!S uv 


- [L] 1 

1 


lyi 
1 1 


1 0 1 
1 1 


! K 


> 

> 


0 T i 
1 


1 1 

iV| - 
1 1 


l-l/Jw'l 
1 1 


! -a^c^w^CL]"’’ 


0 


[H] - w2[Mi]! 


1 1 
IP! 


I [G]P! 



Then, it can be 


transform into: 






l[Kuu] - w2[M] 
1 


- [L] I 

1 


lyi 

1 1 _ 


! — K . V 1 

1 1 


! -a^c^w^CL]''' 


CCH] - w2[Mi3 - [G])| 


1 1 
1 p 1 
1 £_i 


i 0 1 



and: 

1/Z - -jw (1C uv^-U + Kvv.V) 



The user provides the excitation frequencies and, using 
specific boundary condition data, defines the degree of freedom 
that is associated with the applied electrical potential. The 
solver computes, using matrix inversion and Gauss condensation, 
the real and Imaginary parts of the displacement and pressure 
fields. Then, it computes the electrical Impedance. Real and 
imaginary parts are provided In two successive tables. Ten 
successive excitation frequencies can be processed In the same 
run, and, as In the previous section, restart Is available. 

Information given In the previous section and regarding damping 
elements; nearfleid or farfleld conditions, boundary surface and 
interface degrees of freedom remain available for this section and 
are not duplicated. Information given in section 2.9 regarding the 
electrical impedance is also valid. 
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2.13 HARMONIC ANALYSIS OF A PIEZOELECTRIC RADIATING STRUCTURE 
USING THE MIXED METHOD (HAR 5) 



This case is a mixing of the previous cases reported In 
sections 2.10 and 2.12. Using the corresponding notations, the 
matrix equation becomes: 



[Kuu] - w2[M] 


K o 


- [L] ! 

1 


1 Ui 
1 1 


! 0 ! 
1 1 


K T 
t O 


o 

o 


0 T : 

1 


lUo! - 
1 1 


1 Fo ! 
1 1 


-a^c^w^EL]''’ 


0 


[H] - w2[M-,]! 


1 1 
! P! 


![G]P! 



The user provides the excitation frequencies and, using 
specific boundary condition data, defines the degree of freedom 
that is associated with the normal displacement to the Interface. 
The solver computes, using matrix Inversion and Gauss 
condensation, the real and Imaginary parts of the d I sp I acemehTt and 
pressure fields. Real and Imaginary parts are provided In two 
successive tables. Ten successive excitation frequencies can be 
processed In the same run, and, as In the previous section, 
restart I s ava I lab le. 

Information given In section 2.11 regarding damping elements, 
nearfleld or farfield conditions, boundary surface and Interface 
degrees of freedom remain available for this section and are not 
dupl Icated. 
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FINITE ELEMENT MODELING 
OF 

P I EZOELECTR I C TRANSDUCERS 
AUGUST 1987 



CHAPTER 3 

DESCRIPTION OF A DATA FILE 



THE FINITE ELEMENT CODE AT I LA HAS BEEN DESIGNED BY THE ACOUSTICS 
LABORATORY OF THE "INSTITUT SUPERIEUR D' ELECTRON I QUE DU NORD" 
(ISEN. LILLE. FRANCE) FOR THE "GROUPE D'ETUDES ET DE RECHERCHE EN 
DETECTION SOUS-MAR INE" (GERDSM. DCAN. FRENCH NAVY. TOULON). THE 
SOFTWARE ENGINEERING IS PROVIDED BY THE SINAPTEC COMPANY (LILLE. 
FRANCE) . 




I 



Chapter 3. Description of a data file. 



3.1 INTRODUCTION 



This section Is placed first In the data preparation part of 
this manual for pedagogical purposes only. In fact. In most cases, 
an automatic data generation can be processed using the MOSAIQUE 
code and a more simple Initial data file. This automatic data 
generation Is described In the chapter 6. When used first, 
MOSAIQUE provides data files which are Identical to the files 
described hereafter. 

In the general case, an ATILA data file Includes five different 
sect Ions: 

- the header, which contains only comments and describes the Job 
under consideration, 

- a set of entries, which allows, among other things, to describe 
the mesh and the properties of the materials, 

- a set of data, which describes the loading and boundary 
conditions, 

- a set of data, which describes the plane wave part of the 
structure In the case of a mixed model, 

- a set of data, which defines the graphic displays. 

All data, except for the angles, must be given using MKS units or 
any coherent unit system which Is deduced from the MKS system. The 
angle values are always given In degrees. 

In this chapter, we always refer to a simple Illustrative 
example and give the corresponding version of each entry or data 
set. This example Is the modal analysis of an elementary 
ax I symmetr I ca I piezoelectric structure, made up of two ceramic 
rings and two thick electrodes, under electrical ant I resonance 
conditions. This structure has a symmetry plane which is normal to 
the axis and the modes which have to be computed are symmetrical 
wljh respect to this plane. Five full eigenvectors are requested. 
The structure. Its material properties and the corresponding mesh 
are described on the following two pages. Other examples are given 
at the end of the manual. 
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Test example description. 



Geometry. 



2 . 6 . 



///////// 
// AU4G / 
///////// 



4 . 



X5105 



X5105 






2 . 



2 . 



Ground. 



///////// 
// AU4G / 
///////// 



4 . 



Open . 



SYMMETRY Dimensions are in centimeters. 

AXIS 



Material properties. 



AU4G (aluminium alloy). 

.Young's modulus: 0.714E+11 Pa 

.Poisson's ratio: 0.344 

.Density: 2780 kg/m^ 



X5105 (piezoelectric ceramic material). 

.Polarization: axial (up for the top ceramic ring and 
_ down for the bottom ceramic ring). 

.Elastic constants (MKS units, constant electric field): 
S11: 0.114E-10 S33: 0.126E-10 

S12: -0.339E-11 S13: -0.410E-11 

S44: . 0.510E-10 

.Piezoelectric constants (MKS units): 
d33: 0.208E-09 d31 : -0.955E-10 

dl5: 0.491E-09 

.Dielectric constants (MKS units, constant strain): 

EP11: 0.667E-08 EP33: 0.687E-08 
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Mesh . 



X 



11 12 13 



9 



10 




SYMMETRY 

AXIS 



Due to simple symmetry arguments, the mesh is ilmlted to one half 
of the structure. Thus, the OY axis Is the intersection of the 
symmetry plane with the plane of the figure. 
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3.2 THE HEADER 



The header Includes comments which describe the Job under 
consideration. The whole header Is automatically copied on the 
first page of every listing which Is related to the job. The first 
line of the header, generally a title, Is automatically repeated 
ahead of each listing page. To define the header, each line begins 
with the symbol *. The number of lines In a header Is unlimited. 



EXAMPLE. 



+ + + + + + + + + + + + + + + + + + + + + + + + + 
+ *... MODAL ANALYSIS . PIEZOELECTRIC STRUCTURE . 

+ * 

+ »... ANT I RESONANCE CASE 

+ * . . . AX I SYMMETR I CAL MODEL 

+ »... 1/2 OF THE STRUCTURE (SYMMETRY PLANE) 

+ * 

+ »... NUMBER OF ELEMENTS: 2 

+ *... NUMBER OF NODES: 13 

+ *... NUMBER OF DOF: 26 

+ » . . . NUMBER OF MASTER DOF : 23 

+ *. . . DOUBLE PRECISION 

+ * 

+ *... DATA FILE WRITTEN BY: 

+ * . . . DATE 
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3.3 THE ENTRIES 



The various entries which are described hereafter are listed In 
a given order, which corresponds to a classical use but Is not 
mandatory. These entries have to be written using a simple free 
format, described below, which allows an easy reading of the file. 
An entry can be associated to given peculiar words, named 
attributes, which have to be written on the entry line, or to 
given numerical values, named parameters, which have to be written 
on the following line. Entries being read In free format, all the 
attributes and parameters have to be provided, even If they are 
dummy or If their numerical value Is zero. 



3.3.1 DESCRIPTION OF THE FREE FORMAT. 
Control characters. 



. [ ] and [,] separate two words on the same line. They have to be 
used to separate an entry from Its following attribute, to 
separate two following attributes or to separate two following 
parameters. 

. [/] and [■=] separate two following lines (carriage return). As 
well as the carriage return, they have to be used to separate an 
entry from Its following parameter. 

. [?] deletes the preceeding characters of the line, 

. [*] when typed at the beginning of a new line, means that this 
line Is a comment. This line Is not printed on a listing. Such a 
comment has not to be confused with the header. 

. [* ♦] allows the writing of a comment In an entry 

line. This comment Is not printed on a listing. 

Real number syntax. 



The following syntaxes, given as examples, are available: 

1.22E+02 . -1.22E2 

-IE-3 . 1.22 

. ~ - 1 . . -1 

Examples. 

The following sequences are equivalent to define the given 
example: one alphanumeric datum (alpha) on the first line, two 
Integers (1,2) on the second line, three real numbers (1.,2.,3.) 
on the third line: 

. alpha 

1 ,2 

1 . ,2. ,3. 



. aipha-1 ,2 
1 . 2. ,3. 



. alpha/1 2-1 . ,2. 3. 
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3.3.2 THE ENTRY DIRECTORY. 



NAME 


COMP. 


ATTRIBUTES OR PARAMETERS 


P 


ANALYS 1 S 


YES 


STATIC, MODAL or HARMONIC 


7 


CLASS 


NO 


AX 1 SYMMETR 1 CAL , PSTRESS or 
PSTRAIN 


7 


NLOAD 


YES 


NLO 


7 


PRECISION 


NO 


SINGLE or DOUBLE 


8 


LCPDDC 


YES 


NLC 


8 


NODES 


YES 


node coordinates 

(end of the set: dummy line) 


9 


NEWAXES 


NO 


CARTESIAN, CYLINDRICAL or 
SPHERICAL 

XO YO ZO OZZ OYY OXX 


10 


SCALE 


NO 


EX EY EZ 


11 


ELEMENTS 


YES 


Rxxxxx MATERx NGx 

Element topology 

(end of the set: dummy line) 


12 


MATERIALS 


YES 


MATERx 

Properties of the material, 
(end of the set: dummy line) 


14 


GEOMETRY 


NO 


NGx 

Geometrical properties. 

(end of the set: dummy line) 


17 


GEOMETRY 
POLAR 1 ZAT . 


NO 


CARTESIAN, CYLINDRICAL or 
SPHERICAL. 

NGx 

Euler angles, coordinates, 
(end of the set: dummy line) 


18 


REDUCTIONS 


YES 


1 N MDF 


24 


FREQUENCIES 


NO 


Frequency values. 


25 


RADIATION 


NO 


MONOPOLAR or DIPOLAR. 


25 


PR 1 NT 1 NG 


NO 


NPR 


25 


END 


YES 




26 



COMP, means that the entry Is compulsory. 



6 




Chapter 3. Description of a data file. 



3.3.3 DESCRIPTION OF THE ENTRIES. 



In this section, the use of the symbol « In front of the entry 
name means that this entry is compulsory and has to be used in 
every Job. The other entries have to be used only in specific jobs 
or are opt Iona I . 



e ANALYS I S STATIC or MODAL or HARMONIC 

This entry defines the kind of analysis which Is requested In the 
Job (see chapter 2, theoretical background) and Is compulsory. 



Test example. + 

+ ANALYSIS MODAL 
+ 



CLASS AX I SYMMETRICAL or PLSTRESS or PLSTRAIN 

This entry specifies, If necessary, that the analysis corresponds 
to an ax Isymmetr leal , a plane stress or a plane strain model. It 
is only needed If 2D elastic elements are used, but then Is 
compulsory. 



Test example. + 

+ CLASS AX I SYMMETR I CAL 
+ 



e NLOAD 
NLO 

This entry Is compulsory and provides NLO which is: 

. the number of loading cases for a static analysis, 

. 4he number of full eigenvectors which are required In the- case 
of a modal analysis, 

. the number of given frequency values In the case of an In-air 
harmonic analysis, 

. twice the number of given frequency values In the case of an 
harmonic analysis of a radiating structure. 



Test example. + 

+ NLOAD - 5 
+ 
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PRECISION SINGLE or DOUBLE 

This entry Is only used If an eigenvalue computation Is run after 
the assembling and specifies the precision of this computation. 
The default answer is single. 



Test example. + 

+ PRECISION DOUBLE 
+ 



« LCPDDC 
NLC 



The CPDDC array contains the node coordinates (X, Y, Z) and the 
numbers of the node degrees of freedom (translations UX, UY, UZ, 
rotations NX, NY, NZ, pressure P and electrical potential V) . 
These data are listed, with one row per node. In the following 
order : 

X, Y, Z, UX (or P), UY, UZ, NX (or V), NY, NZ 

The LCPDDC entry is compulsory and provides NLC, which is the 
number of terms needed in one row. Thus, In the case of an 
ax isymmetr lea I analysis of an elastic structure NLC equals 5, 
while in an ax isymmetr leal analysis of a piezoelectric structure 
NLC equals 7 and in an analysis of a plate NLC equals 8. if 
useless Intermediate columns are Introduced, due to the given 
value of NLC, the corresponding degrees of freedom will have to 
be deleted later. This case happens, as an example, for the 
degrees of freedom Included In columns 3, 4 and 5 In the plate 
case. 



Test example. -i- 

+ LCPDDC 

+ 7 

+ 
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« NODES 
XI Y1 Z1 
X2 Y2 Z2 
X3 Y3 Z3 



XN YN ZN 
.DUMMY. . 

This entry provides the node coordinates, listed with one row per 
node, In the node number ascending order. The list has to be ended 
by a dummy line and the entry Is compulsory. If, for any reason, 
the node coordinate list is spilt in several successive sets, each 
set requires a new entry NODES and has to be finished by a . new 
dummy line. 



WARNING 



When the model Includes both structure 
and fluid domains, nodes have to be set 
In the following order: 

SOLID - FLUID - RADIATING. 

For ax Isymmetr leal models, the OX axis 
Is necessarily the symmetry axis. 



Test example. 



+ 






+ 


NODES 




+ 


0.00 


0.02 


+ 


0.00 


0.05 


+ 


0.00 


0.08 


+ 


0.01 


0.02 


+ 


0.01 


0.08 


+ 


0.02 


0.02 


+ 


0.02 


0.05 


+ 


0.02 


0.08 


+ 


0.04 


0.02 


+ 


0.04 


0.08 


+ 


0.06 


0.02 


+ 


0.06 


0.05 


+ 


0.06 


0.08 


4 - 






+ 






+ 
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NEWAXES CARTESIAN or CYLINDRICAL or SPHERICAL 
XO YO ZO OZZ OYY OXX 

The initial coordinate system In which the node coordinates are 
given (default system) is cartesian and its origin and axes are 
the absolute origin and axes. The entry NEWAXES allows it to be 
mod I f I ed . 

If the new system Is cylindrical (spherical), the node coordinates 
X, Y, Z, associated to the following NODES entry(ies) become the 
cylindrical (spherical) coordinates R, THETA, Z (R, THETA, PHI). 
Moreover, XO, YO, ZO are the coordinates of the origin of the new 
system, while OZZ, OYY, OXX are the angle values of the Euler 
rotations which transform the Initial axes in the new axes. The 
first rotation OZZ Is around the OZ axis of the Initial system. 
The second rotation OYY Is around the OY axis provided by the 
first rotation. The third rotation OXX is around the OX axis 
provided by the second rotation. Angle values have to be given in 
degrees. 

This entry can be used several times during the description of the 
mesh. Its parameters being always defined with respect to the 
absolute system. To come back to the absolute system. It has to be 
used with the CARTESIAN attribute, the parameters being equal' to 
zero. 



Test example. + 

+ NODES 

+ 0.00 0.02 

+ 0.00 0.05 

+ 0.00 0.08 

+ 0.01 0.02 

+ 0.01 0.08 

+ 0.02 0.02 

+ 0.02 0.05 

+ 0.02 0.08 

+ 

+ NEWAXES CARTESIAN 

+ 0.02 0.0 0.0 0.0 0.0 0.0 . 

+ NODES 

+ 0.02 0.02 

+ 0.02 0.08 

+ 0.04 0.02 

+ 0.04 0.05 

+ 0.04 0.08 

+ 

+ 

+ 



Although this data file is correct, it Is quite evident that It 
has no specific usefulness with respect to this example. 
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SCALE 
EX EY EZ 

if necessary, this entry defines scaie factors (EX, EY, EZ) which 
are appiied to the node coordinates inciuded in the foi iowing 
NODES entry(ies), during the description of the mesh. To switch 
off this scaiing, the SCALE entry has to be used with the 
parameters 0, 0, 0 or 1, 1, 1. 



Test example. + 

+ NODES 

+ 0.00 0.02 
+0.00 0.05 

+0.00 0.08 
+ 0.01 0.02 
+0.01 0.08 
+ 0.02 0.02 
+0.02 0.05 

+ 0.02 0.08 



+ NEWAXES CARTES i AN 

+ 0.02 0.0 0.0 0.0 0.0 0.0 . 

+ SCALE ...; 

+ 2.0 1.0 1.0 

+ NODES 

+ 0.01 0.02 

+ 0.01 0.08 

+ 0.02 0.02 

+ 0.02 0.05 

+ 0.02 0.08 

+ 

+ 

+ 



This exampie is a simpie combination of the entries NEWAXES and 
SCALE, it aiiows the modification of the thickness of the 
eiectrode by simpiy modifying the scaie factor (2.00 in the above 
fiie). 



data file. 
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e ELEMENTS 

RAAAAA MATERA NGA 

N11 N12 N13 .... N1P 

N21 N22 N23 .... N2P 



NK1 NK2 NK3 . . . . NKP 



DUMMY 

RBBBBB MATERB NGB 

N11 N12 N13 N1Q 

N21 N22 N23 N2Q 



NL1 NL2 NL3 NLQ 

DUMMY 



DUMMY 

RZZZZZ MATERZ NGZ 
Nil N12 N13 . . N1R 

N21 N22 N23 . . N2R 

NM1 NM2 NM3 . NMR 

DUMMY 

DUMMY 



This entry provides Information about the elements. It Is 
compulsory. Its parameters are defined set by set, each set 
corresponding to successive elements of the same type. The 
various types of elements, their topologies and their 
corresponding material and geometrical properties are described in 
chapter 4. 

For each set of parameters, the first line contains the name of 
the element type, corresponding to RAAAAA, RBBBBB, ..., RZZZZZ In 
the above entry example, the name of the corresponding material 
data set, denoted MATERx (see the entry MATERIALS for more 
details), and the number of the corresponding geometrical property 
set, denoted NGx (see the entry GEOMETRY for more details). Then, 
each following line contains the actual numbers of the element 
nodes, given In a compulsory order which corresponds to the 
element topology (one or two lines per element, depending upon the 
number of nodes belonging to It, the first line being terminated 
by the character & If a second line Is needed). Finally, the set 
has to be closed by a dummy line. Moreover, the whole entry 
parameter list has also to be closed by a dummy line. Thus, the 
two last lines of this entry are always two dummy lines. 

From the computational point of view, elements are assembled in 
the sequential order of this entry parameter list. Thus, In most 
of the cases, a great deal of attention has to be given to this 
order, which has an Important effect on the array sizes and the 
CPU time requested by the Job. 
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Test example. 



+ 

+ ELEMENTS 

+ R37521 X51 1 

+1 3682457 ... 

+ 

+ R36210 AU4G 2 

+ 6 8 11 13 7 9 10 12 

+ 

+ 



+ 

+ 



The above numberings of the nodes can be easily obtained using the 
corresponding element descriptions In chapter 4. ''X51'' and '‘AU4G'' 
are parameters which refer to the MATERIALS entry. ”1" and ”2" are 
parameters which refer to the GEOMETRY entry. Their meaning is 
discussed In the two following sections. 
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« MATERIALS 
MATER 1 
XII X12 X13 
MATER 2 

X21 X22 X23 . . X2P 



MATERN 

XN11 XN12 XN13 XN14 XN15 XN16 & 

& 

XNM1 XNM2 XNM3 XNM4 XNM5 XNM6 & 
DUMMY LINE 



This entry provides information about the material properties. It 
Is compulsory and has to be closed by a dummy line. Its parameters 
are defined by set. Each parameter set corresponds to a given 
material and is Identified by a name, selected by the user and 
containing less than 8 characters ( MATER 1 , MATER2, .... MATERN in 
the above example). This name links a set of elements from the 
previous ELEMENTS entry to the given material. The list of the 
requested properties is provided In the element description. In 
chapter 4. Thus, as examples, this list contains Young's modulus, 
Poisson's ratio and density in the case of a three-dimensional 
elastic element, and It contains the whole elastic, piezoelectric 
and dielectric tensors In the case of a piezoelectric element. If 
a line is Insufficient to describe the properties of a given 
material. It can be extended using the character &. The order in 
which the physical properties have to be written Is compulsory and 
is given below. 

. elast tc material : 

E NU RO 

Where E Is Young's modulus, NU Is Poisson's ratio and RO Is the 

density. 



. piezoelectric material : 



0.0 


0.0 


RO 


0.0 


0.0 


0.0 


S11 


SI 2 


S13 


S14 


S15 


S16 


S21 


S22 


S23 


S24 


S25 


S26 


S31 


S32 


S33 


S34 


S35 


S36 


S41 


S42 


S43 


S44 


S45 


S46 


S51 


S52 


S53 


S54 


S55 


S56 


S61 


S62 


S63 


S64 


S65 


S66 


dll 


d12 


d13 


d14 


d15 


die 


d21 


d22 


d23 


d24 


d25 


d26 


d31 


d32 


d33 


d34 


d35 


d36 


EP11 


EP12 


EP13 


0.0 


0.0 


0.0 


EP21 


EP22 


EP23 


0.0 


0.0 


0.0 


EP31 


EP32 


EP33 


0.0 


0.0 


0.0 
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where [S] Is the constant electric field elastic, tensor, [d] Is 
the piezoelectric tensor and [EP] Is the constant strain 
dielectric tensor. 

. fluid material : 

COMP .0 RO 

where COMP Is the modulus of compression and RO is the density. 



The material property values of classical or more frequently used 
materials can be written in a special file, named MATER. STD. If It 
is the case, these materials need not be described in the 
MATERIALS entry and their names can be used directly as parameters 
In the ELEMENTS entry. If one of the names belonging to the 
MATER. STD file Is used In the MATERIALS entry, the new property 
values which are listed as parameters will automatically replace 
the old values In the file and will be used by the code during the 
run. In this case, a warning will be printed In the listing file. 
Finally, If a new material name Is defined In the MATERIALS entry, 
this name and the corresponding property values wl IT be 
automatically added to the MATER. STD file. 

The MATER. STD file can be edited by the user for any modification, 
such as adding a new material or correcting old property values. 
Its first line contains a control Integer (generally 80, which Is 
the number of columns), the date of the last edition and Its total 
number of lines. The following lines contain the material property 
values, listed as defined above (1 line for elastic or fluid 
materials, 13 lines for piezoelectric materials). Each of these 
following lines starts with two special fields, containing 
respectively the material name and the number of lines set apart 
for this material (as examples: STEEL 1 or CERAMIC 13 ...). If 
this file is edited, the user has to take care of the final total 
number of lines and. If needed, to correct Its value on the first 
line of the file. An opening error can take place during the 
following run of AT I LA if this number is wrong. 
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Test example. 



+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 



MATERIALS 

AU4G 

• 714E+11 .344 2780 

X51 

0.0 0.0 7350. 0.0 0.0 0.0 

0.114E-10 -0.339E-11 -0.410E-11 
-0.339E-11 0.114E-10 -0.410E-11 

-0.410E-11 -0.410E-11 0.126E-10 

0.0 0.0 0.0 0.510E-10 0.0 0.0 & 
0.0 0.0 0.0 0.0 0.510E-10 0.0 & 
0.0 0.0 0.0 0.0 0.0 0.296E-10 & 
0.0 0.0 0.0 0.0-0.491E-09 0.0 & 
0.0 0.0 0.0 0.491 E-09 0.0 0.0 & 
-0.955E-10 -0.955E-10 0.208E-09 

0.667E-08 0.0 0.0 0.0 0.0 0.0 & 
0.0 0.667E-08 0.0 0.0 0.0 0.0 8. 
0.0 0.0 0.687E-08 0.0 0.0 0.0 . 



0.0 0.0 0.0 & . 

0.0 0.0 0.0 & . 

0.0 0.0 0.0 & . 



0.0 0.0 0.0 & . 



+ 

+ 
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GEOMETRY 

NG1 

XII X12 
NG2 

X21 X22 . . X2P 



NGM 

XM1 XM2 . . . XMQ 

DUMMY LINE 



This entry provides information about the geometrical properties. 
It must be closed by a dummy Hne. Its parameters are defined set 
by set, each parameter set corresponding to a given geometry and 
being Identified by a number selected by the user (NG1 , 
NG2,...,NGM In the above example). For a given geometry, the list 
of the requested properties Is provided In the element 
description. In chapter 4 . Thus, as examples, this list contains 
the thickness In the case of a plate element, and It contains the 
radius of curvature In the case of a damping element. The order In 
which the geometrical properties have to be written Is compulsory. 



Test example. 



In this case, the correct entry for geometrical properties Is the 
GEOMETRY POLARIZATION entry which Is described In the next 
section. In fact. Its format Is Identical to the GEOMETRY entry 
format, and thus the example can be considered as common to the 
two entr les. 
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GEOMETRY POLARIZATION STRING 
NG1 

XII X12 
NG2 

X21 X22 . . X2P 



NGM 

XM1 XM2 . . . XMQ 

DUMMY LINE 



This entry provides Information about the geometrical properties 
when piezoelectric elements are present In the mesh, and It Is 
compulsory In this case. The GEOMETRY entry features are all 
relevant to this new entry and must be used for non-plezoelectr Ic 
elements. However, for piezoelectric elements, specific 
parameters have to be added to define the polarization type 
(cartesian, cylindrical or spherical) and or lentat lon.The 
polarization type Is specified by replacing STRING by CARTESIAN, 
CYLINDRICAL or SPHERICAL In the entry name. 



. If STRING-CARTESIAN, the polarization Is homogeneous In the 
whole element and the natural axes 0X-| , OX 2 and OX 3 have to be 

specified with respect to the global axes OX, OY and OZ (OX 3 Is 

the polarization axis In the classical case of a ceramic 
material). To do this, the user has to provide three Euler angles, 
named ALPHA, BETA and GAMMA, which transform OXYZ Into OX 3 X 1 X 2 and 
are defined as follows (see figure 1 ): 

- ALPHA Is the angle of the first rotation, around OZ, which 

transforms the OXYZ system Into a OX'Y'Z system such that OX' 
belongs to the OZX 3 plane, 

- BETA Is the angle of the second rotation, around OY', which 
transforms the OX'Y'Z system Into a OX^Y'Z', 

- GAMMA Is the angle of the third rotation, around OX 3 , which 

transforms the 0 X 3 Y'Z' system Into the OX 1 X 2 X 3 system. Its value 
Is dummy In the particular case of a ceramic material ( 6 mm 
symmetry class). 

The entry parameters are ALPHA, BETA, GAMMA, listed In this 
compulsory order. Attention has to be given to the exact meaning 
of this "cartesian" polarization. In fact. If the polarization Is 
homogeneous In a three-dimensional element. It Is physically 
homogeneous In the corresponding part of the material. On the 
other hand. If the polarization Is homogeneous and Its direction 
orthogonal to the symmetry axis In the case of an ax Isymmetr leal 
element, this polarization Is physically cylindrical In the 
corresponding part of the material. 
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.If STRING-CYLINDRICAL, the polarization Is radial around a given 
axis O'Z' In the whole element. The user has first to specify a 

O'X'Y'Z' system which contains the O'Z' axis. To do this, he 

provides the 0' coordinates (XO, YO, ZO) as well as the Euler 
angles (ALPHA, BETA, GAMMA) which transform the OXYZ system Into 
the O'X'Y'Z' system and are defined as follows (see figure 2a): 

- ALPHA Is the angle of the first rotation, around OZ, which 

transforms the OXYZ system Into a OxyZ system such that Ox 
belongs to the OZX' plane, 

- BETA Is the angle of the second rotation, around Oy, which 
transforms the OxyZ system Into a OX'yz, 

- GAMMA Is the angle of the third rotation, around OX', which 

transforms the OX'yz system into the OX'Y'Z' system. 

it Is self evident that the choice of the O'X' (or O'Y') axis Is 
arbitrary and has to be the simplest. With these new- reference 
axes, the angle of the rotation around O'Z' which brings the O'X' 
axis to coincidence with the natural MX 3 axis, for each point M 
which Is concerned, is directly computed by the code (see figure 
2b). Finally, the user must also provide the angle DELTA which 
brings, for every point M, the MZ' axis Into coincidence with the 
MX 2 natural axis (see figure 2c). The value of DELTA is dummy In 
the case of a ceramic material ( 6 mm symmetry class). The entry 
parameters are ALPHA, BETA, GAMMA, XO, YO, ZO, DELTA, listed In 
this compulsory order. Attention has to be given to the exact 
meaning of this "cylindrical" polarization. In fact. If the 
polarization is cylindrical In a three-dimensional element. It Is 
physically cylindrical In the corresponding part of the material. 
On the other hand. If the polarization is cylindrical around an 
axis normal to the plane of an ax Isymmetr I ca I element and 
Intersecting the symmetry axis at the point 0, this polarization 
Is physically spherical around 0 In the corresponding part of the 
mater lai . 



. If STRING-SPHERICAL, the polarization Is radial around a given 
point O' In the whole element. The user has first to specify a 
O'X'Y'Z' system centered on the origin O'. To do this, he provides 
the 0' coordinates (XO, YO, ZO) and the Euler angles which 
transform the OXYZ system into the O'X'Y'Z' system (ALPHA, BETA, 
GAMMA) and are defined as for the cylindrical case (see figure 
3a). It Is self evident that the choice of these new axes Is 
arljitrary and has to be the simplest. With these new reference 
axes, the angles of the rotations which bring the O'X' axis Into 
coincidence with the natural MX 3 axis, for each point M which is 
concerned, are computed by the code (see figure 2 b). Finally, the 
user has also to provide the angle DELTA which brings, for every 
point M, the MZ' axis Into coincidence with the natural MX 2 axis 
(see figure 2c). The value of DELTA Is dummy In the case of a 
ceramic material ( 6 mm symmetry class). The entry parameters are 
ALPHA, BETA, GAMMA, XO, YO, ZO, DELTA, listed In this compulsory 
order . 
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Test example. 



+ 

+ GEOMETRY POLARIZATION CARTESIAN 

+ 1 

+ 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

+ 2 

+ 1 

+ 

+ 

+ 



As can be seen after "1”, the three Euler angles are equal to- 
zero, since the polarization direction and the global OX axis are 
identical. Moreover, the value of the thickness after "2" is 
dummy, since the model is ax isymmetr lea I . 
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FIGURE 1. "CARTESIAN" POLARIZATION. 



a/ definition of the Euler angles. 





b/ examples. 

Ax isymmetr leal case. Three-dimensional case. 
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I" 

FIGURE 2. "CYLINDRICAL" POLARIZATION. 

a/ definition of the geometry. 

i ^ 



\ 

\ 








b/ definition of the rotation processed by the code. 




c/ definition of the DELTA angle. 
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FIGURE 3. "SPHERICAL" POLARIZATION. 




b/ definition of the rotations processed by the code. 




c/ definition of the DELTA angle. 
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REDUCTIONS 
I N MDF 

This entry defines the frontal reduction when master degrees of 
freedom are requested (modal and harmonic analysis). In this case, 
it Is compulsory. 

. I Is a control Integer. If 1-0, backsubst I tut Ion Is processed to 
provide NLO full eigenvectors from the eigenmode number N. If 1-1, 
backsubst I tut Ion Is processed to provide NLO full eigenvectors the 
e igenfrequencies of which are above N (In this case, N Is an 
Integer value of a threshold frequency, given In Hz). The control 
Integer N Is mainly devoted to the elimination of the rigid body 
modes, which are provided by the eigenmode computation, from the 
backsubst I tut Ion. In some cases. It can be used to restrict the 
backsubst I tut Ion to a specific given mode. Attention has to be 
paid to the fact that only the twenty first eigenvectors can be 
backsubst I tuted , Including the rigid body modes. 

. MDF Is the chosen number of master degrees of freedom which are 
automatically selected by the code. Attention has to be paid to 
the fact that electrical potential degrees of freedom are always 
condensed, except the applied potential In the case of an harmonic 
analysis which Is a master degree forced by the code. Moreover, 
manually selected degrees of freedom can be defined In the set of 
data which specifies the boundary conditions. Finally, the degrees 
of freedom corresponding to damping elements are all master 
degrees and are forced by the code. Then, the sum of the total 
master degree number IDF (automatic masters + applied potential 
degree + degrees from damping elements) and the automatically 
reduced electrical degree number must be lower than the total 
degree of freedom number IDT. 



Test example. + 

+ REDUCTIONS 

+01 23 

+ 

+ 

Thus, all the available degrees of freedom are master degrees and, 
due to the fact that no rigid body mode exists In this case.^ the 
five first full eigenvectors will be provided. 
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FREQUENCIES 

FI F2 F3 FN 

This entry Is only requested for an harmonic analysis but then Is 
compulsory. Its parameters are the values of the frequency for 
which a computation Is needed. The maximum number of these values 
Is 20 If there Is no damping (In-air harmonic analysis) and 10 If 
damping Is taken Into account (In-water harmonic analysis). The 
parameter line can be extended using the character &. 



Test example. 

There Is no test example for this entry since this entry Is only 
concerned with an harmonic analysis. 



RADIATION MONOPOLAR or DIPOLAR 

This entry Is only requested for the modelling of a radiating 
structure but then Is compulsory. 



Test example. 

There Is no test example for this entry since this entry Is only 
concerned with an harmonic analysis of a radiating structure. 



PRINTING 

NPR 

This entry defines the printing level NPR, which must equal to 1, 
2 or 3, and corresponds to the amount of Information available on 
the result listing. If NPR-1, only the computation results are 
printed. If NPR-2, the more Important data are recalled and 
printed together with the computation results. This level Is the 
normal (default) level. NPR-3 corresponds to the trace level and 
Is only used to detect faults. 



Test example. 

For this example, the default case Is the best choice. 
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9 END 

This entry marks the end of the entry list and Is compulsory. 



+ 

+ END 
+ 

+ 



Test example. 
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3.4 THE LOADING AND BOUNDARY CONDITION DATA 



The data which are associated with the ioading and boundary 
conditions and which are described hereafter have to be written in 
fixed formats. Work is in progress to provide these conditions 
using specific entries in the next ATILA version, in aii the data 
i ines, a given degree of freedom is marked by the number D, with 
the foiiowing correspondence: 1 for UX (or P), 2 for UY, 3 for UZ, 
4 for NX (or V), 5 for NY and 6 for NZ. Moreover, lines which are 
paraiiei to a giobai coordinate axis and planes which are normai 
to a giobai coordinate axis are a iso marked by a number , denoted 
P. Then, the correspondence is: 1 for a piane normai to the OX 
axis, 2 for a piane normai to the OY axis, 3 for a piane normai to 
the OZ axis, 4 for a i ine paraiiei to the OX axis, 5 for a i ine 
paraiiei to the OY axis and 6 for a i ine paraiiei to the OZ axis. 



3.4.1 DESCRIPTION OF THE LOADING CONDITION DATA. 



in the case of a static analysis, the lumped applied forces (or 
moments) have to be described using one of the two following 
procedures. The procedure Is selected using a control number, 
named CONT, which has to be written first. In FI 0.3 format: 

CONT 

Then, If CONT-7777.0, the nodes where the forces (or moments) are 
applied are defined by their coordinates and, for each applied 

force (or moment), the following line has to be written (3F10.3, 
215, F10.3): 

X Y Z D..NL0 F 

X, Y, Z are the coordinates of the fulcrum, D Is the application 
direction (1 to 3 for an applied force, 4 to 6 for an applied 
moment, following the above mentioned correspondence rules), NLO 
Is the loading case number (from 1 to the value defined previously 
bythe NLOAD entry), F Is the force (or moment) value. 

If CONT-8888.0, the nodes where the forces (or moments) are 
applied are defined by their number and, for each applied force 
(or moment), the following line has to be Issued (315, F10.3): 

.NODE D. .NLO F 

NODE is the application node number, D, NLO and F have the same 
meanings as In the previous case. 

In each case, the set of loading condition data lines has to be 
closed by a dummy line. 
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WARNING 



if the model Is ax Isymmetr lea I » the force 
(or moment) actual value has to be divided 
by 2»PI. 

In the case of a distributed loading, the 
lumped nodal forces (or moments) have to be 
computed using the corresponding element 
shape functions. 



Test example. 

The test example Is a modal analysis and does not require any 
loading condition. 



3.4.2 DESCRIPTION OF THE BOUNDARY CONDITION DATA. 



The boundary conditions are generally associated with 
mechanical clampings, symmetry planes, electrodes, pressure 
released surfaces... They are defined by specific data lines, the 
classical format of which Is (315, 4X, A1 ) . The available cases 
are arranged In two groups. The first group contains only boundary 
conditions which refer to the global coordinate axes. The second 
group contains boundary conditions which require specific local 
coordinates. 



Boundary conditions defined In global axes. 

* The degree of freedom D Is deleted at the node N. 

N D 

* The degree of freedom D Is deleted for all the nodes which 
belong to a specific plane or line P containing the node N. 

...-N D P 

Thjs condition can be restricted to the solid or the fluid part of 
a mesh In the case of coupled problems by specifying S or F In the 
fourth field and using one of the two following lines: 

...-N D P S 

. . .-N D P F 

* The degree of freedom D Is deleted for all the nodes In the 
mesh. 



D 
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* Degrees of freedom of type D are Identical for all the nodes 
which belong to a specific line or plane P containing the node N. 

...-N D...-P 

* The degree of freedom D at node N Is a master degree of freedom. 

* « »«N« .« 

* Degrees of freedom of type D are master degrees of freedom for 
ail the nodes which belong to a specific line or plane P 
containing the node N. 

...-N...-D P 

* Degrees of freedom of type D are master degrees of freedom and 
Identical for all the nodes which belong to a specific line or 
plane P containing the node N. 

-N -D -P 

* in an harmonic analysis of a piezoelectric structure, the degree 
of freedom associated with the applied potential has to be a 
master degree of freedom and specifically marked, since" Its 
numerical processing requires a special storage. This selection is 
realised using the following data line: 

where N Is a node which belongs to the related electrode. 
Moreover, this marked degree of freedom Is the same for the whole 
electrode. Thus, If this electrode Is a specific line or plane P, 
the data line can be written: 

... ~*N ... ~4 ... ~P ....1%. 

* In an harmonic analysis of a piezoelectric structure using the 
mixed plane wave-finite element method, the degree of freedom 
associated with the displacement D normal to the interface between 
the two domains has to be a master degree of freedom and 
specifically marked, since its numerical processing requires a 
special storage. This selection is realised using the following 
data line: 

....N... — D ......^..1................................... 

where N is the node which corresponds to D. Moreover, this marked 
degree of freedom Is the same for the whole Interface. Thus, If 
this Interface Is a specific line or plane P, the data line can be 
written: 



.-N...-D 



.-P I 
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* The degree of freedom of type D at node N Is Identical to the 
degree of freedom of the same type D at node M. 

N D M 



Boundary conditions defined In local axes. 

Local axes can be requested In two different cases. Firstly, 
the displacement which Is constrained can be In a direction which 
is not parallel to the global axis (figure 1). Secondly, a 
boundary condition has to be applied on a line which Is not 
parallel to a global axis or on a plane which Is not perpendicular 
to a global axis (figure 2). In fact, these two cases are 
generally merged (figure 3). 



Y 



Y 



P 



0 




P 





X 



X 



Figure 1 



Figure 2 



Y 



D 




X 



Figure 3. 
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* To define local axes at a given node N, the user has to Issue 
the following lines (formats 215 and 6F10.3). 



.All 



.A12. 



.A13. 



.A21 



.A22 / 

. . .A23 



All, A12, A13 are the direction cosines of the first local axis 
and A21 , A22, A23 those of the second local axis. Then, degrees of 
freedom at node N are automatically defined by the code In the new 
local axes and they can be constrained using the same data lines 
as for the first group. If the same global axes have to be defined 
for all the nodes which belong to the same line or plane P 
containing the node N, the preceeding data set has to be modified 
by simply substituting for the first line the fol lowing line: 



.-N. . .10. 



» To constrain nodes belonging to a plane which Is not 
perpendicular to a global axis, the corresponding line of the 
first group has to be modified by substituting 9001 to P and 
adding Immediately after a second line which contain^ the 
direction cosines A1 , A2 and A3 of the normal to this new plane 
(format 3F10.3). Thus, If the degrees of freedom of type D of all 
the nodes which belong to this plane have to be deleted, the two 
data I ines are: 



. . .-N .D.9001 



A1 A2 A3 

To constrain nodes belonging to a line which is not parallel to a 
global axis, the corresponding line of the first group has to be 
modified by substituting 9004 to P and adding immediately after a 
second line which contains the direction cosines A1 , A2 and A3 of 
this new line (format 3F10.3). The preceeding example can be 
immedlatiy transposed. 

• Finally, the two preceeding cases can be merged simply by mixing 
the corresponding data lines. Thus, If the preceeding local axes 
have to be defined for all the nodes belonging to the preceeding 
pl^ne, the following data lines have to be Issued: 

. . .-N. . .10.9001 



All 


A12 


A13 


.A1 


A2 


A3 



A21 A22 / 

. . .A23 
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General cooments. 

. This data line set has to be closed by a dummy line. 

. Constraints on different degrees of freedom can be merged on the 
same data line by placing side by side the corresponding Integers 
In the definition of the D Integer. Thus, If the same constraint 
has to be applied to UX, UY and NZ, D has to be equal to 126. 

. Null flux condition Is the natural condition of the finite 
element formulation used In ATILA and Is Implicit. This Is the 
case for the fluid velocity In fluid domains and for the electric 
field In piezoelectric domains. 



Test example. + 

+ 3 

+ . . .-1 . . .14 5 

+ . . .-6 4 5 

+ 

+ 

+ 

The first line deletes the UZ component of the displacement field 
(compulsory, due to the fact than LCPDDC-7). The second line Is 
the symmetry condition along the OY axis. The third line Is the 
open electrode condition on the upper side of the ceramic ring. 
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3.5 THE MIXED MODEL SPECIFIC DATA 



When a mixed finite eiement-piane wave modei is used, the 
finite eiement part is described simpiy foi lowing the guide ilnes 
given in the other sections. This section is devoted to the 
description of the piane wave part of the structure, which is 
assumed to be organised as shown below. 



.///////////. 
.///////////. 
.///////////. 
./ FINITE //. 
./ ELEMENT /. 
./ DOMAIN //. 
.///////////. 


1 

1 

1 • 

1 

1 

1 

1 • • 

1 

1 

1 

1 

1 

1 


PART 1 


1 

• 1 • • 

1 

• 1 

1 

• 1 • • 

1 

1 . PART 3 . 

1 

1 

• 1 


.///////////. 

.///////////. 


1 

1 

1 

1 

1 

1 

1 

PI 


PART 2 


1 

• 1 
1 

• 1 
1 

1 

P2 



The plane wave part Is limited on one side by Its Interface PI 
with the finite eiement domain. This Interface has to be a plane 
and its normal displacement Is assumed to be uniform, it Is 
limited on the other side by another piane face P3, the normal 
displacement of which is free or clamped. Finally, It can be 
composed of one, two or three parts. Part 1 is the active part and 
its presence is compulsory. Part 3 is passive and in series with 
part 1. it can be missing. Part 2 is also passive and in parallel 
with part 1. It can be missing, but its presence requires the 
presence of part 3. In the classical case of Tonpliz transducers, 
part 1 Is the ceramic stack, part 2 the prestress rod and part 3 
the taiimass, but other schemes can be described if they follow 
the previous general organisation. The description of this 
structure is made in the compulsory order: part 1, part 2, part 3. 
Ail the constitutive elements are assumed to be of cylindrical or 
conical shape. 

The data ilnes are the followings. The first line contains an 
index, named IND (format: 15). 

. . IND 
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if IND-1, the structure Is ax I symmetr I ca I . If IND-2, the structure 
is three-dimensional with one symmetry plane, if IND-3, the 
structure Is three-dimensional with two symmetry planes. The 
following lines describe the successive slices of part 1, 2 and 3. 
Their format Is (2I5, 6E10.3). They have to be written under the 
form: 

. . ITO. . ITH TH 0R1 IR1 OR2 / 

... IR2 PHI 



where: 

ITO - NM, where N and M are two Integers 
defined hereafter, 

ITH - 0 or 180 and allows to reverse the 
ceramic polarization direction. Its value Is 
dummy for other materials, 

TH is the thickness of the slice, 

0R1 is the outside radius of the slice first 
face, 

IR1 Is the Inside radius of the slice first 
face, 

0R2 is the outside radius of the slice second 
face, 

IR2 Is the Inside radius of the slice second 
face, 

PHI is the applied potential (beware of the 
sign) , If any. 

N defines the type of the slice: 

N - 1 for an elastic cylinder (part 1 or 3), 

- 2 for a piezoelectric ceramic cylinder 
(part 1), 

- 3 for an elastic cone (part 1 or 3), 

- -1 , -2 or -3 If this slice Is the last 
slice In part 3, 

- 4 for an elastic cylinder (part 2), 

- 5 for a conical cylinder (part 2), 

- -4 or -5 If this si Ice Is the last slice In 
part 2, 

skits 

M defines the slice material : 

M - 1 for 25CD4 steel or X5105 ceramic, 

- 2 for AU4G or X9 ceramic, 

- 3 for brass (UZ39PB2) or X31 ceramic, 

-4 for "light” steel (RO/100) or X37 
ceramic. 

To close this set, the user has to write down the following data 
I Ine (215) : 

. . . .0. . IP3 



34 



Chapter 3. Description of a data file. 



if IP3 Is equal to zero, the normal displacement to P3 Is free, 
while if IP3 Is equal to 1 this normal displacement Is clamped. 



WARNING: With this version, only 39 successive slices can be 
accomodated by the code. 



*. Cylindrical elements can be replaced by elements with 
different cross section shapes (squares, rectangles ...), If the 
plane wave assumption is still correct in these cases. Then, the 
inside and outside radius of the corresponding cylindrical slices, 
as defined In page 34, have to be computed by the user to ensure 
that the cross section area values are the actual values. 

Material property values are provided by two -.DATA 
instructions in a subroutine named TAC. The user can modify them 
by editing the corresponding FORTRAN file. For elastic materials, 
values are listed In the DATA/SM/ Instruction, with two values for 
each material (density and Young's modulus), successive materials 
being put In the M ascending order. For piezoelectric materials, 
values are listed in the DATA/SP/ Instruction, with eleven values 
for each material (density, constant electric field elastic 
constants S 33 , $ 13 , S-|-|, S-| 2 » S 44 , piezoelectric constants d 33 , 
^31* ^15 constant strain dielectric constants EP 33 AND EP-j-j). 
Successive materials are also put In the M ascending order. 
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3.6 THE GRAPHIC DISPLAY SPECIFIC DATA 



With the ATILA code, graphic displays can be obtained using 
different devices, mainly plotters and workstations. At the 
present time, results have been obtained with BENSON plotting 
tables, VERSATEC Xerox tables, Tektronix 4014, 4114, 4115, as weii 
as MASSCOMP and MiCROVAX computers. When plotters have to be used, 
the ATILA run generates a file which contains the required 
information for the plotter and this plotter is generally 
activated in a following job. In this case, graphic display 
specific data have to be added to the data file and are described 
hereafter. When a workstation Is used, the graphic part of the Job 
is interactive and no data which are related to the graphic 
display have to be put in the data file. Nevertheless, during the 
run of the interactive job, questions are submitted to the user 
which correspond strictly to the following Information. 

Graphic displays of the structure in its rest position or in 
deformed shapes can be easily obtained using the following ^t of 
data. For each graphic display, two data lines are requested with 
the respective formats (4F10.3, 315) and (A40, -F10.3) . These lines 
are : 

ALPHA BETA XLTH BRTH. . .NP. .NE1 . .NE2 

TITLE SCALE 



* ALPHA and BETA are angles (degrees) which define the projection 
of the structure on the drawing plane. The structure being 
initially observed by the user in the OY direction, out of this 
user, the OX axis is rotated of the ALPHA angle around the OZ 
axis. Then, the initial OZ axis is rotated of the BETA angle 
around the original OX axis. Finally, the structure Is projected 
in the initial OXZ plane. ALPHA and BETA are positive if the axis 
are moving back of the initial OXZ plane. 

* XLTH and BRTH (in meters) define the space occupied by the 
graphic display, if BRTH equals zero and XLTH Is negative, ^-XLTH 
Is the scale of the display with respect to the actual size of the 
structure . 

* NP defines the type of display, with the following 

correspondences: 

NP - 0 . Exploded view of the mesh. 

NP - 2 . Exploded view of the mesh, with the 
node numbering. 

NP —2 . Exploded view of the mesh, with the 
e i ement number i ng . 
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NP - 3 . Dashed view of the mesh, with the 
degrees of freedom numbering. 

NP - 7 . Simple display of the mesh (full 
I Ines) . 

NP -12 . Display of the mesh with full outside 
lines and 'dashed Inside lines. 

NP -31 to 49 . Graphic display of the deformed 
shape of the structure. In the loading case, 
full eigenvector or excitation frequency order. 

The deformed shape Is In full line while the 
rest position Is in dashed line, if damping is 
taken Into account In the computation, two 
graphic displays have to be Issued for each 
frequency (real and Imaginary part of the 
displacement field). 

For NP - 0, 2, -2, 3, 7 and 12 the display is provided by a 
preprocessor. If NP - 31 to 49, the display Is provided by a 
postprocessor. Several displays can be simultaneously requested In 
each of these two categories. 

* NE1 and NE2 permit getting the graphic display. of a part of the 
structure only. In this case, NE1 is the number of the first drawn 
element and NE2 the number of the last. If these Integers are 
equal to zero, the whole stucture Is displayed. Care must be given 
to the case of fluid domain, since In this case the finite element 
variable Is the pressure which cannot be represented as a 
displacement. So, fluid elements have to be excluded from graphic 
displays. 

* A specific TITLE can be given to each graphic display. This 
title cannot contain more than 40 letters. 

* SCALE Is a multiplying factor for the displacement amplitude. It 
allows adjustment of this amplitude value for sake of clarity. In 
the case of an eigenmode computation, for which the sign of the 
whole displacement field Is not absolutely defined, it allows also 
to reversal of this sign for the sake of comparison between 
results of different computations ( In the case of piece part 
analysis, for example). 
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CHAPTER 4 
ELEMENT LIBRARY 
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Chapter 4. Element library. 



4.1 INTRODUCTION 



This chapter describes the finite elements which are available 
in the ATILA code . For each element, it provides its present 
name, Its name in the previous ATILA version (in parenthesis). Its 
definition, the list of the degrees of freedom which are attached 
to its nodes (translations, rotations, eiectricai potential, 
pressure), its entry parameters (topology, material properties, 
geometrical properties). 

Most of the elements in the ATILA library are Isoparametric. 
Thus, complex structures with curved sides or faces can be easily 
modelled. Nevertheless, best results are obtained when these 
elements have reasonable shapes. Wild distortions lead to 
inaccuracy because the assumptions made within the cod^ for 
strains become unresilstic. Extreme distortion may even cause 
program failure. The user should guard against excessive 
distortion by adhering to- the following guide lines for 2D 
elements or for faces of 3D elements: 

- the angles between adjacent sides of quadrilaterals should be 
between 45 and 135 degrees, 

- the angles between adjacent sides of triangles should be between 
30 and 100 degrees, 

- in elements with curved sides, the radius of curvature of each 
side should exceed the length of the longest side. 



Warning: all data must be given using MKS units or In any coherent 
unit system which Is deduced from the MKS system. 



*. Design and test of several other elements are in progress In 
the i.S.E.N. Acoustics Laboratory. Further Information can be 
provided upon request. 
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4.2 AT I LA ELEMENT DIRECTORY 



NATURE 


GEOM 


TYPE 


N 


NAME 


COMMENT 


P 


ELASTIC 


, HEXA 


3D 


20 


HEXA20E 




4 




PRIS 


3D 


15 


PRIS15E 




5 




QUAD 


2D&AX 1 


8 


QUAD08E 




6 




QUAD 


PLATE 


8 


PLAT08E 




7 




QUAD 


SHELL 


8 


FACE08E 


FACET 


8 




TRIA 


2D&AX 1 


6 


TR 1 A06E 




9 




TRIA 


PLATE 


6 


PLAT06E 




10 




LINE 


SHELL 


3 


SHEL03E 


AXI 


11 






AXI 


5 


TRAN05E 


TRANS 


12 




LINE 


SPRING 


2 


SPRI02E 




13 


PIEZO. 


HEXA 


3D 


20 


HEXA20P 




14 




PRIS 


3D 


15 


PRIS15P 




15 




QUAD 


AXI 


8 


AXIS08P 




16 




QUAD 


2D 


8 


QUAD08P 


PL. ST 


17 




TRIA 


AXI 


6' 


AX 1 S06P 




18 




TRIA 


2D 


6 


TRIA06P 


PL. ST 


19 


FLUID 


HEXA 


3D 


20 


HEXA20F 




20 




PRIS 


3D 


15 


PRIS15F 




21 




QUAD 


AXI 


8 


AXIS08F 




22 




TRIA 


AXI 


6 


AXIS06F 




23 


SPECIAL 


QUAD 


3D 


16 


QUAD16I 


INTER 


24 




TRIA 


3D 


12 


TRIA12I 


INTER 


25 




LINE 


AXI 


6 


AXIS06I 


INTER 


26 




QUAD 


3D 


8 


QUAD08M 


MONOP 


27 




TRIA 


3D 


6 


TR 1 A06M 


MONOP 


28 




LINE 


AXI 


3 


AXIS03M 


MONOP 


29 




LINE 


AXI 


3 


AXIS03D 


DIPOL 


30 



GEOM-.geometr leal shape 
of the element. 

N : number of nodes. 

NAME: name of the element, 
as required for the 
entry. 

P :page number. 



HEXA :hexaedra. 

PR IS : triangular prism. 
QUAD :Quadr I lateral . 

TRIA -.triangle. 

INTER: Interface element. 
MONOP:monopolar damper. 
DIPOL -.dipolar damper. 
PL.ST:plane strain. 
TRANS:transl t Ion element. 
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Chapter 4. Element library. 



4.3 ELEMENT DESCRIPTION 

Descriptions of a variety of elements are provided on the 
fol lowing pages. 
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Chapter 4. Element || 



HEXA20E (R37110). 



DESCRIPTION. 

.HEXA20E Is a twenty node isoparametric hexaedra, used for 
homogeneous isotropic eiastic materia is. 

.The degrees of freedom are, for each node: Uj{, Uy, Uj, 

(3 trans i at ions) . 

TOPOLOGY. 



7 . . . .20 



8 



18 . 19 . 

15 16 

. . . *1 .... 

. 3. ...12. ...4 



Z 13 . 14 . 

. Y .10 .11 

0...X 1....9....2 



ENTRIES. 

."ELEMENTS" entry parameters. 

HEXA20E MATER 

N1 N2 N3 N4 N5 N6 N7 N8 N9 N10 Nil N12 N13 N14 N15 N16 N17 N18 N19 
N20 

where N1,...N20 are the actual node numbers corresponding to the 
integers 1....20 on the above figure. 

."MATERiALS" entry parameters. 

E NU RO 

where E is Young's moduius (Pa), NU is Poisson's ratio and RO is 
the density (kg/m^). 

."GEOMETRY" entry parameters. 

No. 

COMMENTS. 

.Generally, the aspect ratio has to be less than 3. 
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Chapter 4. Element library. 



PRIS15E (R37210). 



DESCRIPTION. 



.PRIS15E Is a fifteen node Isoparametric triangular prism, used 
for homogeneous Isotropic elastic material. 



.The degrees of freedom are, for each node: u^, Uy, Uj 

(3 translations). 



TOPOLOGY. 



6 



15 14 

12 



Z 

. Y 



0 



. X 




3 



10 




1 ... 7 ... 2 



ENTRIES. 

."ELEMENTS" entry parameters. 

PRIS15E MATER 

N1 N2 N3 N4 N5 N6 N7 N8 N9 N10 Nil N12 N13 N14 N15 

Where N1,...N16 are the actual node numbers corresponding to the 

Integers 1,...15 on the above figure. 

."-MATERIALS" entry parameters. 

E NU RO 

Where E Is Young's modulus (Pa), NU Is Poisson's ratio and RO Is 
the density (kg/m^). 

."GEOMETRY" entry parameters. 

No. 



COMMENTS. 



.Generally, the aspect ratio has to be less than 3. 
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i 



QUAD08E (R36210). 



DESCRIPTION. 

.QUAD08E is an eight node Isoparametric quadrilateral, used for 
homogeneous Isotropic elastic material and for plane stress, plane 
strain or ax Isymmetr leal analysis. 

.The degrees of freedom are, for each, node: u^, Uy 

(2 translations). 

TOPOLOGY. 




ENTRIES. 

."ELEMENTS" entry parameters. 

QUAD08E MATER NG 

N1 N2 N3 N4 N5 N6 N7 N8 

where N1,...N8 are the actual node numbers corresponding to the 
integers 1,...8 on the above figure. 

."MATERIALS" entry parameters. 

E NU RO 

where E Is Young's modulus (Pa), NU Is Poisson's ratio and RO is 
the density (kg/m^). 

."GEOMETRY" entry parameter. 

T 

where T Is the thickness of the element (Ignored In the case of 
plane strain and ax Isymmetr lea I analyses). 

COMMENTS. 



.Generally, the aspect ratio has to be less than 3. 
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PLAT08E (R38110). 



DESCRIPTION. 

•PLAT08E Is an eight node flat quadrilateral plate element, used 
for homogeneous Isotropic elastic material. It relies on the 
classical Love-KI rchhof f hypothesis. 

.The degrees of freedom are, for each node: u^, n^j, Oy 
(1 translation, 2 rotations). 

TOPOLOGY. 



Y 



0 . 



. . X 



3 8 .... 4 



6 7 



1 .... 5 .... 2 



ENTRIES. 



."ELEMENTS" entry parameters. 

PLAT08E MATER NG 

N1 N2 N3 N4 NS N6 N7 N8 

where N1,...N8 are the actual node numbers corresponding to the 
Integers 1,...8 on the above figure. 

."MATERIALS" entry parameters. 

E NU RO 

where E Is Young's modulus (Pa), NU Is Poisson's ratio and RO Is 
the density (kg/m^). 

."GEOMETRY" entry parameter. 

T 

where T Is the thickness of the element. 

COMMENTS. 



.Generally, the aspect ratio has to be less than 3. 

.The longitudinal and lateral dimensions of the whole plate have 
to be at least ten times the thickness. 

.This element can be superimposed to the QUAD08E element (plane 
stress condition) to generate a flat shell element (facet). 
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FACE08E (R 44210). || 



DESCRIPTION. 

.FACE08E Is an eight node quadrilateral shell element, used for 
homogeneous Isotropic elastic material. It relies on the classical 
Love-KI rchhoff hypothesis and has to be flat (facet). 

.The degrees of freedom are, for each node: u^, Uy, U2. n^, ny 
(3 translations, 2 rotations). 

TOPOLOGY. 



3 .... 8 .... 4 



Y 



0 . 



. . X 



6 



7 



1 



5 .... 2 



ENTRIES. 



."ELEMENTS" entry parameters. 

FACE08E MATER NG 

N1 N2 N3 N4 N5 N6 N7 N8 

where N1,...N8 are the actual node numbers corresponding to the 
Integers 1,...8 on the above figure. 

."MATERIALS" entry parameters. 

E NU RO 

where E Is Young's modulus (Pa), NU Is Poisson's ratio and RO Is 
the dens I ty (kg/m^) . 

."GEOMETRY" entry parameter. 

T 

where T Is the thickness of the element . 

COMMENTS . 



.Generally, the aspect ratio has to be less than 3. 

.The longitudinal and lateral dimensions of the whole shell have 
to be at least ten times the thickness. 

.This element Is equivalent to the superposition of the QUAD08E 
element (plane stress condition) and PLAT08E element. 
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TRIA06E (R36110). 



DESCRIPTION. 

.TRIA06E is a six node Isoparametric triangle, used for 
homogeneous isotropic elastic material and for plane stress, plane 
strain or ax isymmetr leal analysis. 

•The degrees of freedom are, for each node: u^, Uy 

(2 translations). 

TOPOLOGY. 



3 



Y 



0 . 



. . X 




ENTRIES. 

.’’ELEMENTS" entry parameters. 

TRIA06E MATER NG 
N1 N2 N3 N4 N5 N6 

where N1,...N6 are the actual node numbers corresponding to the 
Integers 1,...6 on the above figure. 

."MATERIALS" entry parameters. 

E NU RO 

where E Is Young's modulus (Pa), NU Is Poisson's ratio and RO Is 
the density (kg/m^). 

."GEOMETRY" entry parameter. 

T 

where T Is the thickness of the element (Ignored In the case of 
plane strain and ax Isymmetr lea I analyses). 

COMMENTS. 



.Generally, the aspect ratio has to be less than 3. 



9 



Chapter 4. Element 



PLAT06E (R38401). 



DESCRIPTION. 

.PLAT06E Is a six node flat triangular plate element, used for 
homogeneous Isotropic elastic material. It relies on the classical 
Love-KI rchhof f hypothesis. 

.The degrees of freedom are, for each node: U 2 , njf, Oy 

(1 translation, 2 rotations). 

TOPOLOGY. 



3 



Y 



0 



. X 




ENTRIES. 

."ELEMENTS" entry parameters. 

PLAT06E MATER NG 
N1 N2 N3 N4 NS N6 

where N1,...N6 are the actual node numbers corresponding to the 
Integers 1,...6 on the above figure. 

."MATERIALS" entry parameters. 

E NU RO 

where E Is Young's modulus (Pa), NU Is Poisson's ratio and RO Is 
the density (kg/m^). 

."GEOMETRY" entry parameter. 

T 

where T Is the thickness of the element. 

COMMENTS. 



.Generally, the aspect ratio has to be less than 3. 

.The longitudinal and lateral dimensions of the whole plate have 
to be at least ten times the thickness. 

.This element can be superimposed to the TRIA06E element (plane 
stress condition) to generate a flat shell element (facet). 
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SHEL03E (R42123). 



DESCRIPTION. 

.SHEL03E Is a three node ax Isymmetr leal thin shell element, the 
formulation of which takes account of Its double curvature. It 
relies on the classical Love-Ki rchhof f hypothesis. 

.The degrees of freedom are, for each node: Uy, n^ 

(2 translations, 1 rotation). 

TOPOLOGY. 



Y 



0 . . . X 



2 .... 3 



ENTRIES. 



. "ELEMENTS''entry parameters. 

SHEL03E MATER NG 
N1 N2 N3 

where N1 , N2, N3 are the actual node numbers corresponding to the 
Integers 1, 2, 3 on the above figure. 

."MATERIALS” entry parameters. 

E NU RO 

where E Is Young's modulus (Pa), NU Is Poisson's ratio and RO Is 
the density (kg/m^). 

."GEOMETRY" entry parameter. 

T R 

where T Is the thickness of the element and R Its radius of 
curvature (around the axis). 



COMMENTS. 



.The longitudinal and radial dimensions of the whole shell have to 
be at least ten times the thickness. 
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TRAN05E (R 43000). 



DESCRIPTION. 

.TRAN05E Is a five node ax isymmetr leal element which ensures the 
transition between quadrilateral or triangular ax Isymmetr I ca I 
elements and the previous thin shell ax I symmetr I ca I element. It 
allows an accurate description of their Interaction. It Includes a 
solid part (nodes 1, 2, and 3) and a shell part (nodes 4 and 5). 

.The degrees of freedom are, for nodes 1, 2 and 3; u^, u« 
(2 translations), and for nodes 4 and 5: Uj(, Uy, n^ (2 
translations, 1 rotation). 



TOPOLOGY. 

Y 



0 . . . X 



ENTRIES. 

. "ELEMENTS"entry parameters. 

TRAN05E MATER NG 
N1 N2 N3 N4 N5 

where N1 , N2 , N3, N4, N5 are the actual node numbers 

corresponding to the Integers 1, 2, 3, 4, 5 on the above figure. 

."MATERIALS" entry parameters. 

E NU RO 

where E Is Young’s modulus (Pa), NU Is Poisson's ratio and RO Is 
the density (kg/m^). 

."GEOMETRY" entry parameter. 

TP 

where T Is the thickness of the element and R Its radius of 
curvature (around the axis). 



COMMENTS. 

.The longitudinal and radial dimensions of the whole shell have to 
be at least ten times the thickness. 
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SPRI02E (R30100). 



DESCRIPTION. 

.SPRI02E Is a two node spring element, which Is only able to 
transmit to Its two limiting nodes a restoring force due to a 
change of Its length. Its classical use Is for the modelling of 
prestress bolts In transducer design. 

•The degrees of freedom are, for each node: Ux» Uy, 

(3 trans lat Ions) . 

TOPOLOGY. 



0 . . . X 




. 2 



ENTRIES. 

."ELEMENTS" entry parameters. 

SPRI02E MATER NG 
N1 N2 

where N1 and N2 are the actual node numbers corresponding to the 
Integers 1 and 2 on the above figure. 

."MATERIALS" entry parameters. 

E .0 RO 

where E Is Young's modulus (Pa) and RO Is the density (kg/m^). 

."GEOMETRY" entry parameter. 

S 

where S Is the cross section area. 

COMMENTS. 

.In axisymmetr leal models, the actual area S has to be divided by 
2’*‘pl. In three-dimensional models. If the rod belongs to N 
symmetry planes. Its actual area S has to be divided by N+1 . 

.This element Is nothing more than a classical spring . Attention 
has to be paid to spurious modes which can be due to the simple 
spring hypothesis. 
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i 

HEXA20P (R37511). 

DESCRIPTION. 

.HEXA20P Is a twenty node Isoparametric hexaedra, designed for any 
piezoelectric material. 

.The degrees of freedom are, for each node: u^, Uy, u^, v 

(3 translations, 1 electrical potential). 

TOPOLOGY. 



7 ... .20 .... 8 



• • * • • 

18 . 19 . 

15 16 

5. ...17. ...6 . 

3. . .12. . .4 



Z 13 . 14 . 

. Y .10 .11 

0...X 1....9....2 



ENTRIES. 

."ELEMENTS" entry parameters. 

HEXA20P MATER NG 

N1 N2 N3 N4 N5 N6 N7 N8 N9 N10 Nil N12 N13 N14 N15 N16 N17 N18 N19 
N20 

where N1,...N20 are the actual node numbers corresponding to the 
Integers 1,...20 on the above figure. 

."MATERIALS" entry parameters. 

The whole e last I c,p lezoe lectr Ic and dielectric tensors (see the 
detailed description of the entry: MATERIALS). 

."GEOMETRY" entry parameters. 

The new origin coordinates and the Euler angles which define the 
natural axes of the material (see the detailed description of the 
entry: GEOMETRY POLARIZATION). 

COMMENTS. 



.Generally, the aspect ratio has to be less than 3. 

.The entry “GEOMETRY" has to be Issued under the special form 
"GEOMETRY POLARIZATION (CARTESIAN, CYLINDRICAL, SPHERICAL)". 
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PRIS16P (R37512). 



DESCRIPTION. 

•PRIS15P Is a fifteen node Isoparametric triangular prism, 
designed for any piezoelectric material. 

.The degrees of freedom are, for each node: u^, Uy, u^, v 
(3 translations, 1 electrical potential). 

TOPOLOGY. 6 



15 



14 



12 



4 



.13 . . 




Z 

. Y 



0 . . . X 



3 



10 



11 



8 



ENTRIES. 

."ELEMENTS" entry parameters. 

PRIS15P MATER NG 

N1 N2 N3 N4 N5 N6 N7 N8 N9 N10 Nil N12 N13 N14 N15 

where N1,...N15 are the actual node numbers corresponding to the 

Integers 1,...15 on the above figure. 

."MATERIALS" entry parameters. 

Th^ whole elast lc,plezoelectr Ic and dielectric tensors (see the 
detailed description of the entry*. MATERIALS). 

."GEOMETRY" entry parameters. 

The new origin coordinates and the Euler angles which define the 
natural axes of the material (see the detailed description of the 
entry: GEOMETRY POLARIZATION). 

COMMENTS. 



.Generally, the aspect ratio has to be less than 3. 

.The entry "GEOMETRY" has to be Issued under the special form 
"GEOMETRY POLARIZATION (CARTESIAN, CYLINDRICAL, SPHERICAL)". 
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AXIS08P (R37521). 



DESCRIPTION. 

.AXIS08P Is an eight node Isoparametric quadrilateral, designed 
for any piezoelectric material and used only for axisymmetr leal 
analysis. 

.The degrees of freedom are, for each node: Uj{, Uy, v 

(2 translations, 1 electrical potential). 

TOPOLOGY. 



3 .... 8 .... 4 

6 7 

Y 

0...X 1....5....2 



ENTRIES. 

."ELEMENTS" entry parameters. 

AXIS08P MATER NG 

N1 N2 N3 N4 N5 N6 N7 N8 

where N1,...N8 are the actual node numbers corresponding to the 
Integers 1,...8 on the above figure. 

."MATERIALS" entry parameters. 

The whole elastic, piezoelectric and dielectric tensors (see the 
detailed description of the entry: MATERIALS). 

."GEOMETRY" entry parameters. 

The new origin coordinates and the Euler angles which define the 
natural axes of the material (see the detailed description of the 
entry: GEOMETRY POLARIZATION). 



COMMENTS. 



.Generally, the aspect ratio has to be less than 3. 

.The entry "GEOMETRY" has to be Issued under the special form 
"GEOMETRY POLARIZATION (CARTESIAN, CYLINDRICAL, SPHERICAL)". 
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QUAD08P (R37501). 



DESCRIPTION. 

.QUAD08P Is an eight node Isoparametric quadrilateral, designed 
for any piezoelectric material and used only for plane strain 
analysis. 

.The degrees of freedom are, for each node: u^, Uy, v 

(2 translations, 1 electrical potential). 

TOPOLOGY. 



3 .... 8 .... 4 
• • 

• • 

• • 

• • 

6 7 

Y 

• • • 

• • • 

• • • 

0...X 1....5....2 



ENTRIES. 

."ELEMENTS" entry parameters. 

QUAD08P MATER NG 

N1 N2 N3 N4 N5 N6 N7 N8 

where N1 ,.. .N8 are the actual node numbers corresponding to the 
Integers 1,...8 on the above figure. 

."MATERIALS" entry parameters. 

The whole elastic, piezoelectric and dielectric tensors (see the 
detailed description of the entry: MATERIALS). 

."GEOMETRY" entry parameters. 

The new origin coordinates and the Euler angles which define the 
natural axes of the material (see the detailed description of the 
entry: GEOMETRY POLARIZATION). 



COMMENTS. 

.Generally, the aspect ratio has to be less than 3. 

.The entry "GEOMETRY" has to be Issued under the special form 
"GEOMETRY POLARIZATION (CARTESIAN, CYLINDRICAL, SPHERICAL)". 
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AXIS06P (R37522). 

DESCRIPTION. 

.AXIS06P Is a six node Isoparametric triangle, designed for any 
piezoelectric material and used only for ax Isymmetr leal analysis. 

.The degrees of freedom are, for each node: u^, Uy, v 

(2 translations, 1 electrical potential). 

TOPOLOGY. 



3 



Y 



6 



5 



0 . 



. . X 




ENTRIES. 

."ELEMENTS" entry parameters. 

AXIS06P MATER NG 
N1 N2 N3 N4 N5 N6 

where N1,...N6 are the actual node numbers corresponding to the 
Integers 1,...6 on the above figure. 

."MATERIALS" entry parameters. 

The whole elastic, piezoelectric and dielectric tensors (see the 
detailed description of the entry: MATERIALS). 

."GEOMETRY" entry parameters. 

The new origin coordinates and the Euler angles which define the 
natural axes of the material (see the detailed description of the 
entry: GEOMETRY POLARIZATION). 



COMMENTS. 



.Generally, the aspect ratio has to be less than 3. 

.The entry "GEOMETRY" has to be Issued under the special form 
"GEOMETRY POLARIZATION (CARTESIAN, CYLINDRICAL, SPHERICAL)". 
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TRIA06P (R 37502). 



DESCRIPTION. 



.TRIA06P Is a six node Isoparametric triangle, designed for any 
piezoelectric material and used only for plane strain analysis. 

.The degrees of freedom are, for each node; Uj{, Uy, v 
(2 translations, 1 electrical potential). 

TOPOLOGY. 



3 



Y 



6 



5 



0 . 



. . X 



1 



4 



2 



ENTR I E$ . 

."ELEMENTS" entry parameters. 

TRIA06P MATER NG 
N1 N2 N3 N4 N5 N6 

where N1,...N6 are the actual node numbers corresponding to the 
Integers 1,...6 on the above figure. 

."MATERIALS" entry parameters. 

The whole elastic, piezoelectric and dielectric tensors (see the 
detailed description of the entry: MATERIALS). 

."GEOMETRY" entry parameters. 

The new origin coordinates and the Euler angles which define the 
natural axes of the material (see the detailed description of the 
entry: GEOMETRY POLARIZATION). 



COMMENTS. 



.Generally, the aspect ratio has to be less than 3. 

.The entry "GEOMETRY" has to be Issued under the special form 
"GEOMETRY POLARIZATION (CARTESIAN, CYLINDRICAL, SPHERICAL)". 
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HEXA20F (R37810). 

DESCRIPTION. 

.HEXA20F is a twenty node Isoparametric hexaedra, used for fluid 
homogeneous media. 

.The degree of freedom Is, for each node: p (the pressure). 
TOPOLOGY. 

7 ... .20 .... 8 



18 



15 



17 



19 



16 



Z 

. Y 



13 



10 



12 . . . 



14 



11 



0 



. X 



1 



9 



2 



ENTRIES. 

."ELEMENTS" entry parameters. 

HEXA20F MATER 

N1 N2 N3 N4 N5 N6 N7 N8 N9 N10 Nil N12 N13 N14 N15 N16 N17 N18 N19 
N20 

where N1....N20 are the actual node numbers corresponding to the 
Integers 1,...20 on the above figure. 

."MATERIALS" entry parameters. 

COMP .0 RO 

where COMP is the modulus of compression (Pa) and RO Is the 
density (kg/m^). 

."GEOMETRY" entry parameters. 

No. 

COMMENTS. 

.Generally, the aspect ratio has to be less than 3. 



I 
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PRIS15F (R37820). 



DESCRIPTION. 

•PRIS15F Is a fifteen node Isoparametric triangular prism, used 
for fluid homogeneous med I a . 

•The degree of freedom Is, for each node: p (the pressure). 
TOPOLOGY. 

6 



Z 

Y 



0 . 



. . X 



15 14 

12 



4» » .13. .5 



10 



11 

9 8 



1 ... 7 ... 2 



ENTRIES. 

.■'ELEMENTS" entry parameters. 

PRIS15F MATER 

N1 N2 N3 N4 N5 N6 N7 N8 N9 N10 Nil N12 N13 N14 N15 

where N1,...N15 are the actual node numbers corresponding to the 

Integers 1,...15 on the above figure. 

."MATERIALS" entry parameters. 

COMP .0 RO 

where COMP Is the modulus of compression (Pa), and RO Is the 
density (kg/m^). 

."GEOMETRY" entry parameters. 

No. 



COMMENTS. 



.Generally, the aspect ratio has to be less than 3. 
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AXIS08F (R37531). 



DESCRIPTION. 

.AXIS08F Is an eight node Isoparametric quadrilateral, used for 
fluid homogeneous media and for ax Isymmetr lea I analysis. 

.The degree of freedom Is. for each node: p (the pressure). 

TOPOLOGY. 



3 .... 8 .... 4 



Y 



0 



. X 



6 7 



1 .... 5 .... 2 



ENTRIES. 

."ELEMENTS" entry parameters. 

AXIS08F MATER 

N1 N2 N3 N4 N5 N6 N7 N8 

where N1....N8 are the actual node numbers corresponding to the 
integers 1....8 on the above figure. 

."MATERIALS" entry parameters. 

COMP .0 RO 

where COMP Is the modulus of compression (Pa) and RO Is the 
density (kg/m^). 

."GEOMETRY" entry parameters. 

No^ 

COMMENTS. 

.Generally, the aspect ratio has to be less than 3. 
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AXIS06F (R37532). 



DESCRIPTION. 

.AXIS06F Is a six node Isoparametric triangle used for fluid 
homogeneous media and for ax Isymmetr lea I analysis. 

.The degree of freedom Is, for each node: p (the pressure). 

TOPOLOGY. 



3 



Y 



0 



. X 



6 



5 



1 



4 



2 



ENTRIES. 

.“ELEMENTS'* entry parameters. 

AXIS06F MATER 
N1 N2 N3 N4 N5 N6 

where N1,...N6 are the actual node numbers corresponding to the 
Integers 1....6 on the above figure. 

.“MATERIALS" entry parameters. 

COMP .0 RO 

where COMP is the modulus of compression (Pa) and RO is the 
density (kg/m^). 

.“GEOMETRY" entry parameters. 

No_. 

COMMENTS . 

.Generally, the aspect ratio has to be less than 3. 
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QUAD16I (R37910). 



DESCRIPTION. 

.QUAD16I Is a sixteen node Isoparametric quadrilateral element 
which has to be used to ensure the matching between solid and 
fluid meshes along their Interface In the case of three- 
dimensional analyses. It Includes 8 solid nodes and 8 fluid nodes 
which have the same coordinates as the solid nodes. It has no 
th I ckness . 

.The degrees of freedom are, for each solid node: u^, Uy, 

(3 translations), and for each fluid node: p (the pressure) . 

TOPOLOGY. 

SOLID MESH 

7 . . . .16 .... 8 
.3. . . .12. . . .4 



14 . 15 . 

Z Y .10 .11 

• • • • • • 

• • mm mm 

.. 5.. ..13.... 6 . 

O...X 1....9....2 



FLUID MESH 



ENTRIES. 

."ELEMENTS" entry parameters. 

QUAD16I 

N1 N2 N3 N4 N5 N6 N7 N8 N9 N10 Nil N12 N13 N14 N15 N16 

where N1,...N16 are the actual node numbers corresponding to the 

Integers 1,...16 on the above figure. 

."MATERIALS" entry parameters. 

No. 

."GEOMETRY" entry parameters. 

No. 

COMMENTS. 



.The nodes 1, 2, 3, 4, 9, 10, 11, 12 are fluid nodes. The nodes 5, 
6, 7, 8, 13, 14, 15, 16 are solid nodes. 

.Attention must be paid to the orientation of the OZ axis, from 
the fluid domain to the solid domain. 
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TRIA12I (R37920). 



DESCRIPTION. 

.TRIM 2 1 Is a twelve node Isoparametric triangular element which 
has to be used to ensure the matching between solid and fluid 
meshes along their Interface In the case of three-dimensional 
analyses. It Includes 6 solid nodes and 6 fluid nodes which have 
the same coordinates as the solid nodes. It has no thickness. 

.The degrees of freedom are, for each solid node: u^, Uy, Uj 
(3 translations), and for each fluid node; p (the pressure) . 

TOPOLOGY. 



SOLID MESH 



6 

.3 



Z 



0 . 



12 

Y . 9 .8 

4 10 ..... 5. 

.X 1 7 2 



FLUID MESH 



ENTRIES. 

."ELEMENTS" entry parameters. 

TRIM 2 1 

N1 N2 N3 N4 NS N6 N7 N8 N9 NIO Nil N12 

where N1,...N12 are the actual node numbers corresponding to the 
Integers 1,...12 on the above figure. 

."MATERIALS" entry parameters. 

No. 



."GEOMETRY" entry parameters. 
No. 

COMMENTS. 



.The nodes 1, 2, 3, 7, 8, 9 are fluid nodes. The nodes 4, 5, 6, 
10, 11, 12 are solid nodes. 

.Attention must be paid to the orientation of the OZ axis, from 
the fluid domain to the solid domain. 
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AXIS06I (R37541). 



DESCRIPTION. 

.AXIS06I is a six node isoparametric eiement which has to be used 
to ensure the matching between soiid and fiuid meshes aiong their 
interface in the case of axisymmetr icai anaiyses. it inciudes 3 
soiid nodes and 3 fiuid nodes which have the same coordinates as 
the soiid nodes, it has no thickness. 

.The degrees of freedom are, for each soiid node: u^, Uy 

(2 trans iat ions) , and for each fiuid node: p (the pressure) . 

TOPOLOGY. 



Y 

SOL i D MESH 



3 .... 6 .... 4 
O...X 1....5....2 



FLUiD MESH 



ENTRIES. 

."ELEMENTS" entry parameters. 

AXIS06I 

N1 N2 N3 N4 N5 N6 

where N1....N6 are the actual node numbers corresponding to the 
Integers 1,...6 on the above figure. 

."MATERIALS" entry parameters. 

No. 



."GEOMETRY" entry parameters. 
No. 



COMMENTS. 



.The nodes 1, 2, 5 
nodes. 

.Attention must be 
the fluid domain to 



are fiuid nodes. The nodes 

paid to the orientation of 
the solid domain. 



2, 3, 6 are soiid 
the OY axis, from 
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QUAD08M (R37930). 



DESCRIPTION. 

.QUAD08M Is an eight node Isoparametric quadrilateral damping 
element used to take account of a radiation condition. It Is 
restricted to the monopolar (spherical) damping condition. It has 
to be spherical and attached to the outside surface of a three- 
dimensional fluid mesh. 

.The degree of freedom is, for each node: p (the pressure). 
TOPOLOGY. 



3 .... 8 .... 4 



Y 



0 



. X 



6 



7 



1 



5 



2 



ENTRIES. 

."ELEMENTS" entry parameters. 

QUAD08M MATER NG 

N1 N2 N3 N4 N5 N6 N7 N8 

where N1,...N8 are the actual node numbers corresponding to the 
Integers 1,...8 on the above figure. 

."MATERIALS" entry parameters. 

COMP .0 RO 

where COMP Is the modulus of compression (Pa) and RO Is the 
density (kg/m^). 

."GEOMETRY" entry parameter. 

R 

where R Is the radius of curvature of the element. 

COMylENTS. 



.Generally, the aspect ratio has to be less than 3. 

.With respect to the orientation of the whole mesh. Its normal has 
to point outwards . 
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TRIA06M (R 37940). 



DESCRIPTION. 

.TRIA06M Is a Six node Isoparametric triangular damping element 
used to take account of a radiation condition. It Is restricted to 
the monopolar (spherical) damping condition. It has to be 
spherical and attached to the outside surface of a three- 
dimensional fluid mesh. 

.The degree of freedom Is, for each node: p (the pressure). 
TOPOLOGY. 



3 



Y 



6 



5 



0 . . 



. X 



1 



4 



2 



ENTRIES. 



."ELEMENTS" entry parameters. 

TRIA06M MATER NG 
N1 N2 N3 N4 N5 N6 

where N1....N6 are the actual node numbers corresponding to the 
Integers 1,...6 on the above figure. 

."MATERIALS entry parameters. 

COMP .0 RO 

where COMP Is the modulus of compression (Pa) and RO Is the 
density (kg/m^). 

."GEOMETRY" entry parameter. 

R 

where R Is the radius of curvature of the element. 

COMMENTS. 

.Generally, the aspect ratio has to be less than 3. 

.With respect to the orientation of the whole mesh. Its normal has 
to point outwards . 
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AXIS03M (R37551). 



DESCRIPTION. 

.AXIS03M Is a three node Isoparametric damping element used to 
take account of a radiation condition. It Is restricted to the 
monopolar (spherical) damping condition. It has to be circular and 
attached to the outside surface of an ax Isymmetr leal fluid mesh. 

.The degree of freedom Is, for each node: p (the pressure). 

TOPOLOGY. 



Y 



0 . . 



. X 



1 .... 3 .... 2 



ENTRIES. 

."ELEMENTS" entry parameters. 

AXIS03M MATER NG 
N1 N2 N3 

where N1 , N2 and N3 are the actual node numbers corresponding to 
the Integers 1, 2 and 3 on the above figure. 

."MATERIALS" entry parameters. 

COMP .0 RO 

where COMP Is the modulus of compression (Pa) and RO Is the 
density (kg/m^). 

."GEOMETRY" entry parameter. 

R 

Where R Is the radius of curvature of the element. 

COMMENTS. 



.To ensure the correct orientation of this element, Its local OX 
axis must be oriented In the positive direction of the global OX 
axis. Thus, the global x coordinate of node 2 must be higher than 
the global x coordinate of node 1 on the above figure. 
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AXIS03D (T37551). 



DESCRIPTION. 

.AXIS03D Is a three node Isoparametric damping element used to 
take account of a radiation condition. It Includes the dipolar 
damping condition. It has to be circular and attached to the 
outside surface of an ax isymmetr leal fluid mesh. 

.The degree of freedom Is, for each node; p (the pressure). 

TOPOLOGY. 




0 . . 



. X 



3 .... 2 



ENTRIES. 

."ELEMENTS"entry parameters. 

AXIS03D MATER NG 
N1 N2 N3 

where N1 , N2 and N3 are the actual node numbers corresponding to 
the Integers 1, 2 and 3 on the above figure. 

."MATERIALS" entry parameters. 

COMP .0 RO 

where COMP Is the modulus of compression (Pa) and RO Is the 
density (kg/m^). 

."GEOMETRY" entry parameter. 

R 

where R Is the radius of curvature of the element. 

COMMENTS . 



.To ensure the correct orientation of this element, its local OX 
axis must be oriented in the positive direction of the global OX 
axis. Thus, the global x coordinate of node 2 must be higher than 
the global x coordinate of node 1 on the above figure. 
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CHAPTER 5 

DESCRIPTION OF AN AT I LA RUN 



THE FINITE ELEMENT CODE AT I LA HAS BEEN DESIGNED BY THE ACOUSTICS 
LABORATORY OF THE "INST I TUT SUPER I EUR D* ELECTRON I QUE DU NORD" 
(ISEN. LILLE. FRANCE) FOR THE "GROUPE D'ETUDES ET DE RECHERCHE EN 
DETECTION SOUS-MARINE" (GERDSM. DCAN. FRENCH NAVY. TOULON). THE 
SOFTWARE ENGINEERING IS PROVIDED BY THE SINAPTEC COMPANY (LILLE. 
FRANCE) . 
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Chapter 5. Description of an AT I LA run. 



5.1 INTRODUCTION 



The general organisation of an ATILA Job Is described in 
section 1.3. This chapter provides the corresponding organisation 
of an ATILA run in the special case of a VAX or MICROVAX computer 
using the VMS operating system. The same organisation Is also 
correct for other computer types but with different procedures. 
Moreover, the part of this chapter which Is devoted to the graphic 
displays Implies that the GKS library Is available and can be 
linked during the corresponding part of the run, the user being 
then working interactively with the help of a graphic terminal 
(Tektronix 4014, VT240 ...). 



1 



Chapter 5» Description of an AT 



5.2 THE MAIN RUN 



To run an ATILA job, the user must have written a data file as 
described in chapter 3. This data file must be named JOBNAME.ATI, 
ATI being the Job type. Then the user must run the automatic 
program generator, named PGEN, using the entry:' 

RUN PGEN 

During this run, several questions are asked by the code which 
allow the specification of the Job name (JOBNAME) and the Job type 
(ATI), the verification of the analysis type which has been 
Identified by the code, the requirement of special data storages 
for a Job restart (only In the case of an harmonic analysis) or a 
postprocessing (graphic displays ...). PGEN provides three files 
named JOBNAME. FOR, JOBNAME.COM and JOBNAME. LST. JOBNAME. FOR 
contains the ATILA main program. Generally, It can be used wrthout 
any modification. Nevertheless, It has to be edited In few special 
cases: the modal analysis of an elastic structure flooded In an 
incompressible fluid (MOD 5), the radiated pressure computation 
for a given displacement field (HAR 3) and the use, by advanced 
users, of new algorithms which are currently under test (MOD 1, 
MOD 4, ...). In these cases, footnotes are available at the end of 
the corresponding sections In chapter 2 which provide all the 
necessary information. JOBNAME.COM Is used to compile and link the 
program. It need not be edited. Finally, JOBNAME. LST contains a 
copy of the data file, as read during the PGEN run, as well as 
Information about the degree of freedom numbering. It allows a 
very efficient check of the data file and mainly of the boundary 
conditions. Generally, the user must edit this list and perform 
the corresponding check. 

After the PGEN run, the user has to compile the main progam and 
link it to the requested libraries. This step Is simply realised 
using the entry: 

«J0BNAME 

which provides the executable (JOBNAME.EXE). Finally, the user has 
to run the program using the entry: 

RUN JOBNAME 

During this run, the code asks for the Job name (JOBNAME) and the 
Job type (ATI) and Information appears on the user's screen which 
describes the current status (number of the element the code is 
currently assembling, name of the solver the code Is currently 
using ...) and provides the timing. Moreover, the code writes down 
several flies named JOBNAME. LST, JOBNAME. WRK, JOBNAME. SY2 and 
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J0BNAME.SY4. JOBNAME.LST contains a partial copy of the data file, 
depending upon the printing level selected by the user In this 
data file, and all the results (elgenfrequencles, displacement or 
pressure fields associated with successive elgenmodes or 
excitation frequencies, electrical Impedances ... ). The 
JOBNAME.WRK and J0BNAME.SY2 files are working arrays the code 
always needs. They contain useful Information at the end of the 
run only If a restart has to be performed. The J0BNAME.SY4 file 
contains several arrays the code needs for postprocessing. 
Particularly, It can provide the PLO array which contains the 
successive displacement and/or pressure fields as well as the 
CPDDC array which contains the node coordinates and the degree of 
freedom numbering. This file can be easily read by advanced users, 
by editing the JOBNAME.FOR file, deleting all the subroutine calls 
and replacing the WRITE entries devoted to the SY4 file by the 
corresponding READ entries. 
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I 



5.3 THE RESTART 



To restart an ATILA Job Is often necessary when the user Is 
performing an harmonic analysis, since more than 10 or 20 
different frequencies are often required In the frequency range of 
interest. The restart allows the computation of the frequency 
response using previously assembled matrices and often saves a 
large amount of CPU time. Obviously, it can only be performed If a 
full run has been realised before and If the required arrays have 
been stored. Then, the process is exactly the same as for a main 
run, except that the user has to select the restart during the run 
of PGEN, when the corresponding question appears on the screen. 



4 



Chapter 5. Description of an ATILA run. 



5.4 PRE AND POSTPROCESSING 



In this version of the ATILA code, the available pre and 
postprocessors are mainly devoted to the graphic display using the 
GKS package, although many specific postprocessors (farfleld 
computation, transmitting voltage response and directivity 
computations and displays, stress computations and displays, 
isovaiue displays ...) are currently under test but are not 
documented for the present time. At the preprocessing level, the 
code can provide graphic displays of the mesh with different types 
of numbering (nodes, elements, degrees of freedom). At the 
postprocessing level. It can provide the same graphic displays as 
well as distorted shapes of the structure which describe the 
resulting displacement fields. Detailed Information about these 
graphic displays are given in section 3.6. To obtain a graphic 
display, the user must make the entry: — 

RUN PGRAPH 

Then, PGRAPH reads the JOBNAME.ATI file and writes down three 
files named DJOBNAME.LST, DJOBNAME.FOR and DJOBNAME.COM. Moreover, 
during the run, questions appear on the screen which are very 
similar to the question which appears during a PGEN run. The 
DJOBNAME.LST file Is Identical to the JOBNAME.LST file which Is 
obtained after a PGEN run and allows a detailed check of the data 
file if PGRAPH Is run at the preprocessing level. The DJOBNAME.FOR 
file contains the main program which Is needed to obtain the 
graphic displays. It need not be edited. Finally, the DJOBNAME.COM 
file Is used to compile and link the program*. Generally, It need 
not be edited. 

After the PGRAPH run, the user must compile the main program 
and link it to the requested libraries. This step is realised 
simply using the entry: 

9DJ0BNAME 

which provides the executable (DJOBNAME.EXE). Finally, the user 
must run the program using the entry: 

RUN DJOBNAME 

During this run, the code asks the Job name (JOBNAME), the Job 
type (ATI) and requires the choice of the pre or postprocessing 
level. The preprocessing level can be used as soon as the 
JOBNAME.ATI file is available. The postprocessing level can only 
be used if a full ATILA run has previously been performed (main 
run or restart) and If the necessary storages have been made. The 
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whole graphic display phase must be performed using a graphic 
terminal and the information about the display type, its 
orientation. Its title ..., are required during the run by 
questions which are written down on the screen. 



I 



i i 

I 



*. The link of the main program to the GKS library is provided 
by the file LNK1D.COM which Is used by DJOBNAME.COM. Before the 
first PGRAPH run, the user must check that this link Is correct, 
taking account of the GKS installation on the corresponding 
computer . 
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5.5 THE AUTOMATIC MESH GENERATION 



The use of the automatic mesh generator MOSAIQUE Is described 
In the following chapter. It requires a JOBNAME.MOS file which Is 
written down by the user. Then the user runs MOSAIQUE with the 
entry: 

RUN MOSAIQUE 

MOSAIQUE provides two files named JOBNAME.LST and JOBNAME.ATI . 
JOBNAME.LST contains information which are only needed If a 
problem appears during the generation. The JOBNAME.ATI file Is the 
data file for the ATILA main run. Depending upon the contents of 
the JOBNAME.MOS file, it can be used Immediately or can be edited 
by the user to add the last entries and data, generally the 
loading or boundary conditions. 
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5.6 SUMMARY 



Automatic mesh 
generat Ion. 



Preprocessing. 




JOBNAME.MOS 



RUN MOSAIQUE 

JOBNAME.LST 

JOBNAME.ATI 





JOBNAME.ATI 

RUN PGRAPH 1 

D JOBNAME.LST <— 
D JOBNAME . FOR *- 
DJOBNAME.COM 4-i 

«DJ0BNAME 1 

D JOBNAME . EXE ^ 

RUN OJOBNAME 



Interactive graphic display 



Main run. 



-p JOBNAME.ATI 
U RUN PGEN 



JOBNAME.LST 

E JOBNAME . FOR 
JOBNAME.COM 

«JOBNAME 

p JOBNAME.EXE 

-U RUN JOBNAME 





JOBNAME.LST *— 
JOBNAME. WRK *— 
JOBNAME. SY2 
JOBNAME. SY4 
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Restart . 




JOBNAME , AT I 
JOBNAME.WRK 
JOBNAME. SY2 
JOBNAME. SY4 

RUN PGEN 

JOBNAME. LST 
JOBNAME. FOR 
JOBNAME.COM 

6J0BNAME 

JOBNAME . EXE 

RUN JOBNAME 




JOBNAME. LST ♦- 
JOBNAME.WRK 
JOBNAME. SY2 — 
JOBNAME. SY4 *- 



Postprocessing. 




JOBNAME. AT I 
JOBNAME. SY4 

RUN PGRAPH ■ 

D JOBNAME. LST 
D JOBNAME. FOR 
DJOBNAME.COM 

«DJOBNAME 

DJOBNAME.EXE 

RUN DJOBNAME 








Interactive graphic display 



Due to their large size, the files .LST, .WRK, .SY2, .SY4 are not 
stored as successive versions If they are provided several times 
with the same Job name but are automatically deleted by the code. 
If, for any reason, the user wishes to save one of the previous 
version before a new run with the same Job name. It has to rename 
the corresponding file. 
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CHAPTER 6 

AUTOMATIC MESH GENERATION 



THE FINITE ELEMENT CODE AT I LA HAS BEEN DESIGNED BY THE ACOUSTICS 
LABORATORY OF THE "INSTITUT SUPERIEUR D* ELECTRON I QUE DU NORD" 
(I SEN. LILLE. FRANCE) FOR THE "GROUPE D' ETUDES ET DE RECHERCHE EN 
DETECTION SOUS-MARINE" (GERDSM. DCAN. FRENCH NAVY. TOULON). THE 
SOFTWARE ENGINEERING IS PROVIDED BY THE SINAPTEC COMPANY (LILLE.* 
FRANCE) . 



Chapter 6. Automatic mesh generation. 



6.1 INTRODUCTION 



As mentioned in section 1.3.2 and in the introduction of the 
third chapter, partial regularity of the mesh can often allow the 
use of an automatic generation, which provides node coordinates 
and element topologies. In the ATILA code, this work Is performed 
by a specialized preprocessor named MOSAIQUE. Then, the user has 
only to define a gross splitting, built with super elements, and 
to select their automatic splitting Into finite elements 
(Figure 1). From the numerical point of view, each super element 
is transformed Into Its natural axes, using the Isoparametric 
transformation. Then, It Is split into finite elements, following 
the corresponding data provided by the user, and these finite 
elements are transformed back to the global axes using the r^erse 
isoparametric transformation. Finally, redundant nodes are 
eliminated by searching nodes which belong to a same given small 
2D or 3D window. 



Entries for MOSAIQUE are provided by a specific Input file 
written by the user and named JOBNAME.MOS. At the end of the 
MOSAIQUE run, the output file which contains the generated node 
coordinates and topologies Is named JOBNAME.ATi (see chapter 5). 
Moreover, entries and data which are required by ATILA but are not 
necessary for the mesh generation can be also read by MOSAIQUE, if 
present in the JOBNAME.MOS file, and be copied In the JOBNAME.ATI 
output file. Thus, starting with a simple input file 
(JOBNAME.MOS), MOSAIQUE is able to provide the whole ATILA data 
file (JOBNAME.ATi). Finally, It has to be noted that, because the 
number of nodes which have to be defined for MOSAIQUE Is often 
small, this procedure allows the user to modify very easily 
geometrical parameters (length, width ...) as well as the mesh 
steps. 



_Most of the entries and data contained In the MOSAIQUE.. data 
file are identical to the ATILA entries and data. However, the 
NODES and ELEMENTS entries are modified, while two new entries can 
be used. These entries and data are described in section 6.2. 
Moreover, the splitting of a super element Into finite elements 
can be obtained using different frames. The frames which are 
retained in MOSAIQUE are described In section 6.3. 
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FIGURE 1. SCHEMATIC DESCRIPTION OF THE GENERATION PROCESS. 
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6.2 MOSAIQUE ENTRIES AND DATA 



6.2.1. THE HEADER. 

Information contained In section 3.2 Is correct for the header 
of a MOSAIQUE data file. This header Is reproduced as the header 
of the resulting JOBNAME.ATI file. Thus, It will be the header for 
the ATILA Job. 



6.2.2. THE ENTRIES. 

The ATILA entries ANALYSIS, CLASS, NLOAD, PRECISION, LCPDDC, 
MATERIALS, GEOMETRY, GEOMETRY POLARIZATION, REDUCTIONS, 
FREQUENCIES, RADIATION and PRINTING, If provided In the MOSAIQUE 
data file, are simply read and copied In the resulting JOBNAME.ATI 
file. Obviously, they are not used by the mesh generation 
procedure. The END entry has to be provided to MOSAIQUE. As well 
as for ATILA, It defines the end of the entry list. It will also 
be copied In the resulting JOBNAME.ATI file. 

The NODES entry Is used to provide the coordinates of the nodes 
which define the super elements. Its syntax Is exactly the same as 
for the ATILA data file and the node coordinates have to be listed 
in the node number ascending order associated to the gross mesh. 
After the MOSAIQUE run, this NODES entry wi i I be copied In the 
JOBNAME.ATI file with the coordinate list which Is associated with 
the generated nodes. The new node numbering will be the node 
numbering provided by this list. It has to be noted that the 
numbering of the nodes belonging to the Initial gross mesh Is 
modified during the generation. Moreover, as well as for ATILA, 
nodes have to be set In the compulsory order: solid, fluid, 

radiating. Finally, the entry NEWAXES and SCALE can be used to 
define the node coordinates for the gross mesh as well as they are 
used for ATILA. 

As described in chapter 3, the ELEMENTS entry Is divided In 
sets, each set being related to a given element type . Successive 
sets have to be separated by a dummy line and the whole entry must 
also be closed by a dummy line. However, with respect to ATILA, 
the content of a set Is modified In MOSAIQUE by replacing each 
topology line of the type: 

Nil N12 N13 . . . .NIP 

by a set of lines of the type: 

Nil N12 N13 . . . . NIP 

-1 ITYP NK NL NM NN I 
1X1 1X2 1X3 ... . IXNK 
IY1 IY2 IY3 . . . lYNL 
IZ1 IZ2 IZ3 IZM4 



3 



Chapter 6. Automatic mesh gen 



where: 

- Nil, N12, N13, .... NIP are the actual numbers of the 

corresponding super element nodes, 

- the second line (-1 ITYP NK NL NM NN I) provides compulsory 
Information about the splitting of this super element, 

- ITYP refers to a given frame (see section 6.3), 

- NK Is the number of generated elements along the OX natural axis 
of the super element, 

- NL Is the number of generated elements along the OY natural axis 
of the super element (If necessary), 

- NM Is the number of generated elements along the OZ natural axis 
of the super element (if necessary), 

- NN Is a dummy field, kept for special uses, 

-I Is a control Integer which must be equal to 0 If the 
splittings are regular and 1 If at least one splitting Is 
irregular , 

- 1X1, 1X2, 1X3, ..., IXNK define the splitting along the OX 
natural axis of the super element. 

- IY1, IY2, IY3, ..., lYNL and IZ1, IZ2, IZ3 IZNM provide 

the same Information, respectively for the OY and OZ natural 
axes of the super element. 



Then, along the OX natural axis of the super element, the 
corresponding side Is divided Into NK different parts, the lengths 
of which are respectively equal to: 

IX1.L/IX, IX2.L/IX, IX3.L/IX IXNK.L/IX 

where L Is the total length, IX Is given by: 

IX - 1X1 + 1X2 + 1X3 + ... + IXNK 

and 1X1, 1X2, 1X3, ...,IXNK have to be Integers. The same 

splitting rule Is applied for the OY and OZ natural directions. If 
I Is equal to 0, the three lines which are devoted to the 
description of the splitting are not needed. If I Is equal to 1, 
these three lines must be provided, even If one or two of them 
correspond to a regular splitting. 

In the super element topology, a mid side node can always be 
omitted If the corresponding side Is straight. Nevertheless, If a 
topology line Is Incomplete, It has to be completed with zeros, at 
least until the last defined node. As an example. If the mld-sIde 
nodes 12, 14 and 21 can be omitted In the topology line 1 3 5 9 12 
14 16 21, this line can be written: 

125900 16 0 or 125900 16 

In opposition to the AT I LA case, when solid and fluid elements are 
used together, the elements must be provided In the order: solid. 
Interface, fluid, damping. This order can be modified by the user 
In the JOBNAME.ATI file after the MOSAIQUE run, if needed. 
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Finally, two new entries are available for MOSAIQUE with 
respect to ATILA. The first Is used to provide Information related 
to the generation of the damping elements and must be written In 
the following form: 

SORTFLUID 
NFLU MAC NSH 

where NFLU Is the actual number of the first node belonging to the 
boundary surface of the gross mesh, NAC Is the actual number of 
the node which is the boundary surface centre (acoustic centre) 
and NSH is a dummy field kept for special uses, if the boundary 
surface centre does not belong to the mesh and cannot be 
identified with a node, NCA has to be equal to zero and the origin 
of the global axes must be this centre. This entry Is compulsory 
If damping elements are used. The second new entry Is used to 
modify the size of the 2D or 3D window which allows the 
elimination of redundant nodes. It has to be used If redundant 
nodes are still present In the JOBNAME.ATI file or If actual nodes 
have been deleted. It has to be written under the form: 

MERGE 

XM _ 

where XM will be the new size of the window (}(M Is a length and 
must be provided in meters). This entry has only to be used If the 
above mentioned problems appear, which could be the case If the 
mesh Is very Irregular. Then the best way of chosing the XM value 
Is a trial and error method. 



In fact, the element type is related simultaneously to the 
super element and the generated finite elements. This situation Is 
simple when the super element and the generated elements 
correspond to the same geometrical shape (splitting of a 
quadrilateral super element Into QUAD08E finite elements, 
splitting of a triangular super element into AXIS06P finite 
elements ...) but can lead to severe difficulties if geometrical 
shapes are different (splitting of a quadrilateral super element 
Into TRIA06E finite elements, as an example), in this last case, 
the user must provide an element type which corresponds to the 
geometrical shape of the super element. Then, MOSAIQUE writes down 
this wrong element type In the JOBNAME.ATI file and the user has 
to modify this type before using the data file. Work is In 
progress to avoid this problem In the next ATILA version. 
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6.2.3. THE DATA. 

Data can be simply added to the entries In the JOBNAME.ATI 
file, after the automatic mesh generation. This procedure Is 
classical, since the node numbers which are necessary to write 
down the data are only available after the generation. However, in 
special cases (multiple runs with the same topology and different 
dimensions, for example), the user may wish to write down these 
data In the JOBNAME.MOS file. Then, the rules are the same as for 
an ATILA data file and the corresponding lines are simply read by 
MOSAIQUE and copied Into the JOBNAME.ATI file. 
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6.3 THE GENERATING FRAMES 



The splitting of a super element In simple finite elements can 
be realised In various ways, depending upon the shape of this 
super element and the type of the requested finite elements 
(quadrilateral 2D elements with 8 nodes, 20 node hexaedras, 6 node 
ax isymmetr Icai interfaces ...). The available types of splitting, 
named frames In the previous sections, are marked by the ITYP 
control integer and described in the following lines. 



ITYP - 1 

A linear super element is split Into 2 node linear finite elements 
(SPRI02E) following the diagram-. 



NK elements (NK « 2) 



The super element Is defined by nodes 1 and 2 and the control 
integers NL and NM are dummy. 



ITYP ■ 2 

A linear super element is split Into 3 node damping elements 
(AXIS03M or AXIS03D) following the diagram: 



1 ... 0 ... 0 ... 0...2 

NK elements (NK - 2) 



The super element is defined at least by the nodes 1 and 2 but 
generally also by its mid node, due to the curved shape of the 
damping boundary. The control integers NL and NM are dummy. 
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■i 

I 

ITYP - 3 

A linear super element Is split Into 6 node interface elements 
(AXIS06I) following the diagram: 



Fluid domain 

0 ... 0 ... 0 ... 0...0 

1 ... 0 ... 0 ... 0...2 

Solid doma I n 

NK elements (NK - 2) 



The super element is defined by the nodes 1 and 2. its mid node is 
required only If the interface is curved. The control integers NL 
and NM are dummy. The orientation of the normal to the Interface 
from the fluid to the solid being compulsory for the gen^ated 
Interface finite elements, the orientation of the super element Is 
also compulsory. Nodes 1 and 2 are defined In the above order If 
the Oz natural axis Is the outward normal to this paper sheet. 
Then, the Ox natural axis is oriented toward the left (attention 
must be devoted to this point, since this numbering rule which 
applies to the super element Is the opposite of the rule which 
applies to the finite elements). 



ITYP - 4 

A quadrilateral super element is split Into 8 node quadrilateral 
elements (QUAD08E, PLAT08E, FACE08E, AXIS08P. QUAD08P. AXIS08F or 
QUAD08M) following the diagram: 



NL 

elements 
(NL - 2) 



NK elements (NK - 3) 
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The super element Is defined by the nodes 1, 2, 3 and 4. Mid side 
nodes are only required If the sides are curved. The control 
integer NM Is dummy. 



ITYP - 5 



A quadrilateral super element Is spilt Into 6 node triangular 
elements (TRIA06E, PLAT06E, AXIS06P. TRIA06P, AXIS06F or TRIA06M) 
following the diagram: 



• « 

0 0 

NL 

elements 0 . . 0 . . 
(NL - 3) . 

• « 

0 0 

• • 

« « 

0 . . . 0 
• • 

• • 

0 0 

* * 

• * 
1 . . 0 . . 

NK 



• ••O... 0 •••4 

• « • • 

• • • • 

0 0 0 0 

• • • . 
• • • • 

0 .. 0 .. 0 .. 0..0 

• • • • 

• • • • 

0 0 0 0 

• • • • 

• • • ♦ 

«** 0 ««« 0 ...o 

• • • • 

• • • • 

0 0 0 0 

• • • * 

. • • • 

0..0..0..0..2 

ilements (NK - 3) 



The super element is defined by the nodes 1, 2, 3 and 4. Mid side 
nodes are only required If the sides are curved. The control 
integer NM is dummy. This case Is related to the comments at the 
end of section 6.2.2 and care must be paid to the element type (as 
an example, the element type can be QUAD08E In the MOSAIQUE data 
file and TRIA06E in the ATILA data file). 



ITYP - 6 

A triangular super element Is split Into 6 node triangular 
elements (TRIA06E, PLAT06E, AXIS06P, TRIA06P, AXIS06F or TRIA06M) 
following the diagram: 
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3 

• • 

• • 

0 0 
• • 

• • 

0 . . 0 . .0 

NL . . . . 

elements 0000 

(NL - 3) . . . 

• • • • 
0 .. 0 .. 0 .. 0..0 
• • • • • • 

• • • • • • 

0 0 0 0 0 0 
• • • • • • 

• * • • • • • 

1 .. 0 .. 0 .. 0 .. 0 .. 0..2 

NK elements (NK - 3) 

The super element is defined by the nodes 1, 2 and 3. Mid side 
nodes are only required If the sides are curved. The control 
integer NM Is dummy. 



ITYP - 9 

A quadrilateral super element is split into 16 node quadr i iaterai 
interface elements (QUAD16I) following the diagram: 




NK elements (NK - 2) 
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The super element Is defined by the nodes 1, 2, 3 and 4. Mid side 
nodes are only required If the sides are curved. The control 
Integer NM Is dummy. Attention must be paid to the orientation of 
this super element, the normal to the Interface being always 
oriented from the fluid toward the solid. 



ITYP ■ 10 



An hexaedric super element Is split Into 20 node hexaedrlc finite 
elements (HEXA20E, HEXA20P or HEXA20F) following the diagram: 



7. .0. .0. .0. .8 



0 * « 0 * « 0 * « 0 . *0 



0 . . 0 . . 0 . . 0 . .0 



0* .O* .O. .O. *0 



1 .. 0 .. 0 .. 0..2 

NK elements (NK - 2) 



m 

elements 
(NM - 3) 



NL 

elements 
(NL - 2) 



The super element Is defined by the nodes 1, 2, 3, 4, 5, 6, 7 and 
8. Mid side nodes are only required If the sides are curved. 
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tTYP - 11 

A prismatic super element with a triangular basis Is spilt Into 15 
nodes prismatic finite elements (PRIS15E, PRIS15P or PRIS15F) ' 
following the diagram: 
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0 



0 

NM 

. elements 
0 (NM - 3) 



4.. 0..0..0..5 0 0 

• • « • • • • 

• • # • • • • 

0 O 0 0 0 0 0 

• • • • • • • 

• • • • • • • 
0..0..0..0..0 0 3 

• • • • • • 

« • • • • • 

0 0 0 0 0 0 

• • • • « • 

NL 

0..0..0..0..0 0 elements 

(NL - 2) 

0 0 0 0 

• • • • 

• • • • 

1 .. 0 .. 0 .. 0..2 



NK elements (NK - 2) 



The super element is defined by nodes 1, 2, 3, 4, 5 and 6. Mid 
side nodes are only required If the sides are curved. 



ITYP - 12 

An hexaedric super element Is spilt Into 15 node prismatic finite 
elements (PRIS15E, PRIS15P or PRIS15F) following the diagram: 
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7 . . 0 



0 . . 8 



5 . . 0 




0 0 



. 0 



0 0 0 0 



1 .. 0 .. 0 .. 0 .. 0..0 
NK elements (NK - 3) 



0 0 0 0 



NM 

elements 
(NM - 2) 



NL 

elements 
(NL - 2) 



0 0 



The super element Is defined by the nodes 1, 2, 3, 4, 5, 6, 7 and 
8. Mid side nodes are only required If the sides are curved. This 
case Is related to the comments at the end of section 6.2.2 and 
care must be paid to the element type (as an example, the element 
type can be HEXA20E In the MOSAIQUE data file and PRIA15E In the 
ATILA data file). 



ITYP - 14 

A linear super element Is split into 3 node ax Isymmetr leal shell 
finite elements (SHEL03E). The Information related to ITYP - 2 Is 
also correct for this case but the numbering of the generated 
nodes and the topology of the generated elements take account of 
the specific numbering given for this shell In section 4.3 
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